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1 . 0  INTRODUCTION 


Propeller-hull  interaction  arises  in  two  related  but  distinct  ways. 
Unsteady  forces  and  moments  developed  by  the  propeller  operating  in  the  spa¬ 
tially  nonuniform  hull  wake  and  transmitted  to  the  hull  through  the  propel¬ 
ler  shafting  are  known  as  bearing  forces.  Periodic  pressure  forces  induced 
on  the  hull  due  to  the  hydrodynamic  action  of  unsteady  propeller  velocity 
and  pres.=ure  field  are  termed  surface  forces.  The  increase  of  shaft  output 
power  and  restrictive  profile  of  current  single-screw  designs  give  rise  to 
vibration  and  noise  problems.  In  the  past,  propeller- induced  ship  hull  vi¬ 
bration  has  been  handled  primarily  as  in-service  engineering  modifications. 
Recently,  however,  many  articles  have  appeared  on  the  topic  of  ship  vibra¬ 
tion  and  noise  arising  from  pressure  excitation  (1),  (2).  A  fund  of  theo¬ 
retical  knowledge  on  the  various  hydrodynamic  aspects  of  the  interaction 
problem  has  emerged  from  the  collective  works  of  groups  of  investigators 
such  as  Breslin  et  al.  (3),  Cox  et  al.  (4),  Kaplan  (5),  Noordzij  (6),  Huse 
(7)  and  others.  In  spite  of  these  advances  in  the  understanding  of  the 
flow,  there  is  no  satisfactory  answer  to  the  question  of  whether  some  par¬ 
ticular  wake,  together  with  a  given  propeller  configuration,  will  cause  vi¬ 
bration  problems  at  a  given  speed. 

Briefly,  lifting  surface  methods  have  been  developed  to  calculate  the 
unsteady  flow  field  around  an  Isolated  propeller  operating  in  the  prescribed 
velocity  field  of  the  hull  wake  (8)  and  (9).  The  introduction  of  the  hull 
into  the  known  flow  field  of  the  propeller  results  in  a  time -dependent 
boundary  condition  on  the  hull  surface  which  must  be  satisfied  by  placing 
time -dependent  singularities  on  the  hull  surface.  The  force  and  moments  on 
the  hull  follow  by  an  application  of  the  extended  Legally  theorem.  Vorus 
(10)  has  suggested  an  alternative  to  the  calculation  of  the  hull  diffraction 
potential  in  the  presence  of  the  propeller  onset  flow.  He  has  shown  that 
the  unsteady  forces  and  moments  on  the  hull  can  be  calculated  by  combining 
the  solution  for  the  potential  about  a  bare  hull  in  simulated  motion  with 
the  solution  for  the  isolated  propeller  flow  field.  In  both  of  these  ap¬ 
proaches  the  feedback  of  the  hull  singularities  on  the  propeller  flow  field 
is  neglected. 

The  mechanism  producing  large  pressure  pulse  excitation  involves  a  com¬ 
plicated  interaction  between  cavitating  flow  over  the  propeller  blades  and 
rapidly  changing  velocities  associated  with  nonuniform  wakes.  The  develop¬ 
ment  of  methods  to  model  the  generation  and  volume  variation  of  a  transient 
cavity  is  an  important  prerequisite  for  calculating  propeller-hull  interac¬ 
tion.  Although  the  methodology  for  calculating  propeller-hull  interaction 
in  an  approximate  manner  is  available,  the  prediction  of  cavity  geometry  is 
not  well  developed. 

The  objective  of  the  present  work  is  to  further  the  understanding  of 
the  propeller-hull  interaction  process  by  numerical  simulation.  Specifical¬ 
ly,  the  aim  was  to  couple  the  three-dimensional  inviscid  flow  procedures  de¬ 
veloped  at  Analytical  Methods,  Inc.  (AMI)  with  a  companion  viscous  flow  pro¬ 
gram  at  Scientific  Research  Associates,  Inc.  During  the  term  of  the  con¬ 
tract  the  previously  proposed  investigation  was  revised  in  the  light  of  dis¬ 
cussions  held  with  the  Office  of  Naval  Research.  The  major  change  proposed 
by  the  Navy  was  an  emphasis  on  the  hub  and  tip  vortex  modeling  for  a  shroud¬ 
ed  propeller  in  nonuniform  flow  without  the  participation  of  the  viscous 
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flow  program  from  Scientific  Research  Associates.  The  redirection  of  ef¬ 
forts  at  AMI  lead  to  a  comprehensive  computational  approach  aimed  at  solving 
specific  areas  of  concern  to  the  Navy.  Special  emphasis  was  given  to  pro¬ 
peller  analysis  in  prescribed  hull  wakes  and  to  the  solution  of  the  complete 
inviscid  hull  interaction  problem.  The  basic  approach  usrd  a  superposition 
of  the  solution  for  several  harmonics  of  the  prescribed  nonuniform  velocity 
distribution.  This  capability  was  developed  in  the  program,  MPROP,  a  marine 
propeller  version  of  VSAERO.  The  case  of  a  propeller  in  nonuniform  flow  was 
also  analyzed  using  an  existing  time -dependent  program,  USAERO,  developed  at 
Analytical  Methods,  Inc.  These  two  approaches  are  particularly  suitable  for 
certain  types  of  nonuniformity  of  flow.  The  frequency-based  approach  pro¬ 
vides  the  basic  features  of  the  Interaction  in  a  simple,  cost  effective  way. 
The  time -dependent  solution  is  computationally  less  efficient,  but  it  may  be 
the  only  way  to  solve  an  arbitrary  inflow  problem;  also  it  has  the  potential 
for  a  complete  simultaneous  analysis  of  the  hull/propeller  mutual  interac¬ 
tion.  Both  these  approaches  have  to  be  applied  to  analyze  the  interaction 
problem. 


2.0  PROPELLER  ANALYSIS 


Predicting  the  performance  of  marine  propellers,  even  in  isolation  or 
in  the  ideal  conditions  found  in  most  laboratory  test  environments,  is  a 
challenging  task.  The  high  levels  of  skew  and  twist,  the  high  solidity  of 
the  multi-blade  configurations  now  being  studied,  and  a  basic  inability  to 
model  leading-edge  flows  (important  for  cavitation  studies)  have  forced  the 
designer  away  from  the  simpler  lifting  line  and  lifting  surface  analyses 
that  have  laid  the  foundations  of  the  theoretical  understanding  of  propeller 
flows.  This  has  led  to  the  introduction  of  first  generation  surface  singu¬ 
larity  models  where  the  full  blade  geometry,  including  thickness,  is  repre¬ 
sented. 

Even  these  methods  were  unable  to  properly  model  the  flows  around  the 
blade  root  sections  although  leading-edge  effects  were  now  more  correctly 
represented.  The  use  of  simple  source  or  vortex- lattice  models  of  the  sec¬ 
tions  has  led  to  problems  associated  with  mass  flow  loss  through  the  surface 
(leakage).  This  was  particularly  problematic  when  the  flows  were  channel¬ 
like,  as  in  the  interblade  passages  at  the  roots  of  the  blades  of  today's 
high  solidity  marine  propellers. 

They  were  also  generally  unable  to  allow  for  the  mutual  deformation  of 
the  blade  wakes  and,  more  importantly,  although  some  limited  representation 
of  the  tip-edge  shedding  was  possible,  they  did  not  have  the  capability  to 
model  the  rolling  up  of  the  edge  wakes  and  their  subsequent  development  over 
the  tip  surface.  Finally,  the  ability  to  include  any  viscous  flow  effects 
was  severely  limited  and  prediction  and  modeling  of  regions  of  separated 
flow  was  not  possible. 

These  problems,  the  poor  representation  of  the  root  channel  flows,  the 
wake  development  and  the  viscous  flow  effects,  have  been  resolved  with  a 
second  generation  panel  method,  program  VSAERO  (11)  (for  Vortex  Separation 
AERCWynamics) .  Developed  originally  for  generalized  aero -/hydrodynamic  con¬ 
figuration  modeling  purposes,  the  program  uses  a  compound  source/doublet 
singularity  to  represent  the  surface  and  applies  a  Dirichlet  boundary  condi¬ 
tion  to  solve  for  the  unknown  potential.  This  contrasts  with  earlier  panel 
methods  which  applied  the  Neumann  boundary  conc^tion  to  solve  for  a  locally 
unknown  velocity.  Program  VSAERO  is  also  able  to  calculate  the  development 
of  the  wake  using  an  iterative  procedure  and  uses  a  streamline  tracking  cal¬ 
culation,  together  with  a  quasi-three-dimensional  boundary  layer  analysis  to 
determine  the  viscous  flow  behavior. 

The  program  differs  from  other  second  generation  panel  methods  (Ref.  12 
is  typical)  in  that  a  low-order  (i.e.,  flat  panel)  representation  of  the 
surface,  with  all  the  benefits  of  ease  of  configuration  assembly,  has  been 
retained.  Other  methods,  while  retaining  a  simple  singularity  model,  have 
gone  to  a  higher-order  surface  representation.  This  has  resulted  in  high 
levels  of  effort  to  prepare  the  input  deck  even  for  the  simple  cylindrical 
surfaces  of  aircraft 

A  special  version  of  program  VSAERO,  named  MPROP,  has  been  developed 
for  the  marine  propeller  application.  The  major  changes  to  the  program  were 
a  move  to  a  helical  wake  relaxation  scheme  and  the  addition  of  the  blade 
rotations  to  the  conventional  translation  motions  in  the  original  program. 
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In  this  regard,  it  has  much  in  common  with  another  special  version  of  VSAERO 
developed  to  model  hovering  rotors  or  aircraft  propellers  in  axial  flight 
(13).  The  major  difference,  of  course,  is  the  geometric  representation  of 
the  blades  with  the  marine  version  developing  the  highly  skewed,  twisted  and 
cambered  sections  found  in  present  practice.  The  modeling  has  been  updated 
to  enforce  the  Kutta  condition,  especially  for  propeller  shapes  without  a 
sharp  trailing  edge. 

All  of  the  earlier  marine  propeller  analysis  methods  made  some  assump¬ 
tions  regarding  the  axisymmetry  of  the  flow  and  here  MPROP  was  no  different. 
Consequently,  application  of  the  analyses  was  limited  to  uniform  onset  flow, 
flows  far  removed  from  the  normal  operating  environment  of  the  marine  pro¬ 
peller  . 

The  following  assumptions  are  made  in  the  analysis  of  bodies  moving 
through  a  mass  of  fluid. 

1.  The  fluid  is  at  rest  when  the  body  is  stationary. 

2.  The  fluid  boundaries  are  sufficiently  far  away  from  the  body  and  their 
movement  is  negligible. 

3.  The  fluid  is  homogeneous  and  is  of  uniform  density. 

4.  The  effects  of  viscosity  are  negligible  everywhere  except  inside  ideal¬ 
ized  thin  layers  across  which  the  velocity  potential  is  discontinuous. 

5.  The  speed  of  motion  of  the  body  is  negligible  compared  to  the  speed  of 
sound  in  the  fluid. 

All  of  the  above  assumptions  are  conventional  and  a  variety  of  problems 
important  to  a  marine  engineer  can  be  adequately  modeled  within  the  frame¬ 
work  of  the  above  assumptions.  Section  5  discusses  how  to  analyze  the  mo¬ 
tion  of  a  body  in  a  fluid  which  is  not  at  rest.  The  modeling  of  the  nearby 
boundaries  is  explained  in  Section  2.4.  Section  3  deals  with  flow  situa¬ 
tions  where  the  effects  of  viscosity  may  not  be  negligible.  The  fifth  con¬ 
dition  above  may  be  relaxed  provided  the  flow  is  subsonic  everywhere.  In 
other  words,  all  steady  flow  problems  arising  in  moderately  subsonic  speeds 
are  included.  No  assumptions  or  approximations  regarding  the  shape  of  the 
body  are  made.  There  are  no  limitations  in  the  analysis  which  preclude  the 
application  of  the  code  to  blunt,  abrupt,  or  uneven  body  geometry.  In  the 
case  of  massive  flow  separation,  the  analysis  must  be  repeated  with  dif¬ 
ferent  parameters  or  with  a  different  code  as  dictated  by  the  nature  of  the 
problem.  Boundary  layer  analysis  of  patches  with  roughness  is  possible. 
Cavitation  control  devices  such  as  mass  injection  or  suction  can  be  modeled. 
The  basic  theory  is  contained  in  two  documents  already  submitted  to  the  Navy 
(13),  (14),  A  brief  summary  is  presented  for  completeness. 


Tne  velocity  vector  in  potential  flow  can  be  expressed  as  the  gradient 
of  a  scalar  potential,  The  mathematical  conditions  to  be  satisfied  by  ^ 

are ; 

(a)  -  0  everywhere  in  the  flow  field  except  inside  the  idealized  wake 

layers , 

(b)  70  •  n  —  on  the  surface  of  the  body.  V  is  nonzero  only  when  mass 

is  injected  or  drawn  in  to  induce  the  §esired  flow  pattern,  and 

(c)  70-0  sufficiently  far  away  from  the  body  for  the  disturbance  to  van¬ 

ish. 

Let  w  denote  the  wake  surface  of  bodies  with  nonzero  circulation  such 
that  there  is  a  discontinuous  jump  in  velocity  potential  across  5t.  In  the 
domain  excluding  the  potential  barrier  of  the  wake,  the  flow  is  non-circula- 
tory  and  the  velocity  potential  at  any  point,  P,  in  the  flow  following  Lamb 
(Ref.  7,  p.  59,  Eq.  (7))  is 

0p(x,y,z,t)  -  -1/Un  JJ  1/r  d4>/3n  ds  +  1/Un  Jj  0  d/dn  (1/r)  ds  (2.1) 

sW 


where 
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y(Xp  -  x^)"  +  (yp  -  y^)2  +  (Zp  -  z^)2 

and  Xj^,  y^^  and  z,  represent  a  point  on  the  body  surface.  The  integration 
must  extend  over  alt  points  on  the  body  and  wake  surface.  0  is  the  distur¬ 
bance  potential  due  to  body  motion  and  vanishes  at  infinity.  When  point  P 
is  on  the  surface,  the  locally  singular  kernel,  1/r,  must  be  replaced  by  ap¬ 
propriate  limiting  value  and  should  be  excluded  from  the  region  of  integra¬ 
tion. 


Discretization  of  S  into  N  panels  assuming  that  0  is  constant  on  each 
of  the  panels  and  30/3n  takes  the  value  of  prescribed  normal  velocity,  pro¬ 
vides  the  following  matrix  equation  for  the  unknown  velocity  potential  on 
the  surface,  0  . 


[A]  [0^]  -  (0^1  (2.2) 

Elements  of  matrix  A  are  geometric  quantities  defined  by  JJ  d/dn  (l/r)ds, 

and  0JJ  is  the  velocity  potential  induced  at  a  control  point  by  the  sum  of 
all  constant  strength  sources  of  known  strength.  The  former  may  be  recog¬ 
nized  as  the  velocity  potential  induced  at  a  control  point  by  a  unit 
strength  doublet  whose  axis  is  normal  to  the  panel  surface.  The  above  ma¬ 
trix  equation  is  the  basic  equation  and  is  sometimes  referred  to  as  an  in¬ 
ternal  Dirichlet  formulation  for  the  velocity  potential. 
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2 . 1  Uniform  Flow 


In  the  case  of  a  propeller,  we  have  a  steady  linear  velocity  of  advance 
along  the  shaft  axis  and  rotation  with  a  constant  angular  velocity,  u,  about 
the  shaft  axis.  The  normal  velocity  at  a  control  point,  P,  due  to  a  com¬ 
bination  of  translation  and  rotation  is  given  by  V  •  n  +  (w  x  7)  •  n,  where 
7  is  the  position  vector  of  the  point,  P.  The  change  in  the  expression  for 
the  source  strength  is  one  major  change  in  the  boundary  condition.  The 
other  change  is  the  fact  that  the  trailing  wake  has  a  component  of  swirl 
velocity  added  to  it. 

In  developing  MPROP  for  analyzing  marine  propellers  several  technical 
problems  arose  which  had  to  be  solved  before  smooth  surface  pressures  could 
be  calculated. 

1.  The  standard  geometry  interpolation  package  led  to  excessive  distortion 
and  displacement  in  the  location  of  panel  edges  along  a  column.  This 
was  corrected  by  adding  a  new  interpolation  routine  which  constructs 
panel  edge  coordinates  based  on  direct  interpolation  of  pitch,  camber 
and  thickness. 

2.  The  existing  wake  models  for  wings  and  helicopter  rotor  blades  were  un¬ 
suitable  for  modeling  marine  propeller  wakes  due  to  major  differences 
in  spanwise  loading.  A  new  wake  model  has  been  synthesized  based  main¬ 
ly  on  Min's  experimental  studies  at  M.I.T.  (16). 

3.  The  severe  correction  to  pitch  due  to  large  skew  or  warp  caused  diffi¬ 
culty  for  the  existing  scheme  of  satisfying  the  Kutta  condition,  es¬ 
pecially  near  the  hub.  An  alternative  formulation  of  the  Kutta  condi¬ 
tion  was  introduced  along  with  changes  in  the  wake  routines  to  restore 
correct  loading  behavior. 

4.  In  the  existing  panel  model  the  singularity  distributions  are  placed  on 
the  mean  plane  of  the  quadrilateral.  On  quadrilaterals  whose  sides  are 
not  coplanar  (i.e.,  skewed  panels),  the  comer  points  are  equidistant 
from  the  constructed  mean  plane  (two  above  and  two  below).  This  model, 
therefore,  lea  s  gaps  between  the  singularity  distributions  on  neigh¬ 
boring  panels;  a;  i  result,  the  influence  coefficient  acting  at  nearby 
control  points  can  be  in  error.  The  most  often  seen  effect  of  this  is 
after  wake  relaxation  when  some  of  the  wake  panels  become  highly  skew¬ 
ed.  If  part  of  the  surface  paneling  is  close  by,  such  as  in  close  vor¬ 
tex/surface  interaction  for  tip  vortices  or  propeller -hub  Juncture  vor¬ 
tices,  then  significant  local  errors  occur  in  the  doublet  solution  and, 
hence,  in  the  local  velocities  and  pressures.  Skewed  surface  panels 
are  also  a  source  of  error,  but  the  associated  error  can  be  alleviated 
by  repaneling.  The  most  obvious  treatment  is  to  divide  the  skewed 
panel  into  two  triangles  having  a  common  edge  in  one  of  the  diagonals 
of  the  quadrilateral.  Figure  2.1.  This  treatment  has  proven  adequate 
for  wake  panels  with  little  cost  penalty. 
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Figure  2.1.  Two  Triangles  in  Version  C.0  Wake  Panels. 


on  the 

sphere 

where 


1.0 


Figure  2.2.  Sphere  in  Non-uniform  Flow. 


5. 


Practical  propeller  designs  allow  thick  and  abrupt  trailing  edges. 
Correct  modeling  of  such  edges  poses  a  technical  problem.  The  usual 
practise  of  relaxing  the  sectional  circulation  until  the  satisfaction 
of  the  Kutta  Condition  at  the  trailing  edge  is  not  enough.  Such  an  ap¬ 
proach  results  in  unacceptable  errors  in  overall  performance  predic¬ 
tions.  To  correct  this  problem  a  Y-wake  attachment  has  been  developed 
which  leads  to  satisfactory  solutions  for  the  cases  analyzed.  The  use 
of  two  parallel  wake  surfaces  leads  to  numerical  errors  in  the  calcu¬ 
lated  value  of  circulation.  In  the  Y-wake  model  the  wakes  from  the  up¬ 
per  and  lower  surfaces  are  brought  to  a  line  and  merged  as  a  single 
sheet.  Such  an  approach  eliminates  the  numerical  errors. 


2 . 2  Nonuniform  Inflow 

MPROP  allows  the  specification  of  a  general  field  of  inflow  velocity. 
A  nonuniform  velocity  field  may  be  specified  in  three  different  ways.  The 
three  components  of  the  velocity  field  may  have  an  arbitrary  radial  varia¬ 
tion,  planar  variation  or  spatial  variation.  Cartesian  polar  coordinates 
may  be  used  to  specify  the  inflow  field.  A  three-dimensional  interpolation 
scheme  has  been  devised  to  obtain  the  velocity  components  at  any  given  point 
from  a  field  specification  of  velocity  components. 

2.2.1  Comparison  with  Exact  Solution 

In  order  to  verify  the  modeling  of  nonuniform  inflow,  a  test  case  was 
constructed.  There  is  a  very  special  exact  nonuniform  flow  about  a  sphere 
for  which  the  Laplacian  field  can  be  constructed  assuming  such  flow  is  pos¬ 
sible  without  viscosity.  The  flow  and  its  inviscid  solution  are  shown  in 
Figure  2.2. 

For  such  a  flow  the  normal  velocity  induced  on  the  sphere  is 

3^/37  -  -U  cos  $ 

-  -  1/2  (3-5  cos*  9)  cos  9 

-  PjC/*) 

where  n  -  cos  9  and  P,(m)  is  a  Legendre  polynomial  of  the  third  order.  By 

invoking  the  property  of  zonal  harmonics  it  is  possible  to  write  down  the 

expression  for  the  velocity  potential  on  the  sphere  in  similar  fashion  to 
solving  for  uniform  flow  past  a  sphere  (6),  (see  Article  9(1  in  Ref.  15).  ^ 

on  the  sphere  -  R/8  (5  cos*  J  -  3f  cos  tf). 

The  solution  was  also  calculated  using  VSAERO.  Figures  2.3,  2.4  and 

2.5  compare  the  velocity  potential,  X  component  of  velocity  and  total 
velocity  of  the  calculation  with  the  analytical  solution. 
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EXACT  SVaOTlON 
VSAERO 
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Figure  2.3.  Comparison  of  Velocity  Potentials  on  a  Sphere. 


EXACT  SOLUTION 


Figure  2.4.  Comparison  of  X  Components  of  Velocity  on  a  Spher 


The  above  calculations  were  perfonned  with  the  nonuniform  inflow  pre¬ 
scribed  on  a  radial  line.  Two  additional  calculations  were  made  on  a  sphere 
with  nonuniformity  prescribed  along  a  two-dimensional  plane  and  a  three-di¬ 
mensional  volume,  respectively.  The  doublet  solution  for  three  different 
nonuniform  inputs  is  presented  in  Figures  2.6(a),  (b) ,  and  (c) .  All  three 
forms  are  thus  verified.  A  fourth  form  of  representing  the  nonuniform  flow 
is  also  provided  in  terms  of  polar  x,  r  and  9  components.  This  option  is 
useful  when  the  geometry  of  the  velocity  scan  is  annular. 
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Figure  5.6.  Concluded 


3.0  HUB  AND  SHROUD  INTERACTION 


MPROP  was  modified  to  accept  an  axisymmetric  duct,  an  axisymmetric  hub 
and  an  arbitrary  number  of  blade  rows.  The  modification  allows  for  some 
components  to  rotate  while  keeping  other  components  stationary.  For  ex¬ 
ample,  the  propeller  and  hub  may  rotate  while  the  duct  and  other  components 
are  held  stationary.  In  the  following  the  special  changes  to  the  program, 
MPROP,  which  help  in  modeling  the  hub  and  shroud  are  discussed. 


3 . 1  Basic  Modeling 

In  applying  MPROP  to  model  close  surface  interactions,  three  major  dif¬ 
ficulties  were  encountered.  These  are  related  to  off-body  velocity  calcula¬ 
tions  wake  relaxation  and  streamline /boundary  layer  analysis. 

Off-body  velocity  calculations  have  had  difficulties  in  the  past  when 
performed  near  a  surface  or  inside  a  duct  or  between  nearby  components.  New 
schemes  have  been  installed  to  improve  this  calculation.  A  velocity  calcu¬ 
lation  is  regarded  as  "safe"  when  beyond  a  panel  "size"  from  the  panel  cen¬ 
ter.  The  initial  size  of  a  panel  is  the  average  distance  between  opposite 
mid-points.  Once  the  surface  velocities  have  been  evaluated,  the  distances 
between  opposite  mid-points  are  now  weighted  by  the  local  velocity  compo¬ 
nents  in  these  directions.  Thus,  large  aspect  ratio  panels  with  flow  pre¬ 
dominantly  across  the  small  width  have  a  reduced  near- field  "size"  relative 
to  the  old  model.  In  the  new  velocity  scan  all  "safe"  velocities  along  a 
scan  line  are  computed  first.  All  points  falling  within  the  near-field  zone 
of  a  panel  are  marked  and  their  velocities  are  derived  later  by  interpola¬ 
tion  through  the  "safe"  set.  The  surface  velocity  values  on  the  panels  in¬ 
tersected  by  the  scan  line  are  included  in  the  interpolation.  A  warning  is 
printed  if  an  isolated  near- field  point  is  detected  with  no  external  "safe" 
velocity  information.  The  interpolations  are  performed  independently  within 
each  flow  region  in  cases  where  a  scan  line  has  several  intersections. 
Points  falling  inside  a  closed  body  are  not  treated  except  that  the  first 
point  inside  a  surface  is  transferred  to  the  surface  and  takes  on  the  local 
panel  value.  At  present,  there  is  a  scheme  installed  to  interpolate  for  the 
velocity  at  the  intersection  point- -the  nearest  panel  value  is  used.  With 
reasonable  panel  density  this  should  not  cause  a  major  discrepancy;  however, 
should  it  appear  necessary,  a  local  interpolation  scheme  has  also  been  in¬ 
stalled  to  improve  this  calculation. 

New  scan  options  have  been  installed  in  MPROP  to  allow  panel-to-panel 
and  point- to-panel  scan  lines.  These  give  good  results  except  very  close  to 
the  surface.  The  boundary  layer  velocity  profile  has  been  incorporated  in 
the  off-body  velocity  scan  to  give  a  realistic  velocity  profile  near  the 
surface . 

In  wake  relaxation  calculations  on  complicated  configurations  Involving 
multiple  wakes,  wake  lines  from  different  wakes  have  been  observed  to  pair 
up  to  form  a  doublet  couple  which  proceeds  to  track  alarmingly  off  to  the 
side.  While  this  usually  does  not  seriously  affect  the  solution  on  the  up¬ 
stream  components,  it  can  change  local  wake  roll-up  and  affect  downstream 
components.  A  scheme  has  been  installed  in  MPROP  which  allows  the  formation 
of  multiple  wake  lines.  If  any  two  wake  lines  approach  within  a  certain 
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distance,  they  are  identified  together.  A  bundle  can  include  a  number  of 
wake  lines  and  new  members  may  join  at  the  downstream  wake  grid  planes. 
Only  the  trajectory  of  the  wake  line  with  the  smallest  subscript  in  each 
bundle  is  computed;  the  movement  vectors  are  then  transferred  to  the  other 
members  of  the  bundle. 

Clearly,  the  action  of  the  new  bundle  scheme  supersedes  the  earlier 
coincident  vortex  line  routine  and  makes  it  unnecessary  to  input  coincident 
vortex  line  information  where  wake  lines  become  coincident  at  some  down¬ 
stream  wake  grid  plane. 

The  small  distance  parameter  used  as  a  criterion  for  bundle  formation 
is  presently  formed  on  the  basis  of  the  segment  lengths  on  the  wake  attach¬ 
ment  line  projected  in  the  y-z  plane.  A  reasonable  size  parameter  has  been 
established  on  the  basis  of  a  number  of  test  cases. 

The  wake  geometry  was  contained  in  251  x  31  arrays.  The  maximum  number 
of  wake  lines  was  251  and  the  maximum  number  of  grid  planes  was  31.  The 
small  number  of  wake  grid  planes  is  inadequate  for  the  relaxation  of  tip 
vortices.  In  MPROP,  the  arrays  were  converted  to  single-subscript  arrays, 
so  that  for  those  cases  which  use  fewer  than  251  wake  lines,  the  number  of 
wake  grid  planes  available  is; 

Maximum  grid  planes  <  7,781/(Total  number  of  wake  lines)  (i) 

Maximum  grid  planes  <  155  (ii) 

A  streamline  procedure  has  been  installed  which  is  more  robust  and  in¬ 
cludes  previously  unavailable  options,  such  as  automated  streamline  starts 
and  laminar/turbulent  starts;  that  is,  boundary  layer  analyses  along 
streamlines  which  do  not  start  at  stagnation  points.  Also,  a  significant 
effort  has  been  expended  in  improving  the  reliability  of  the  boundary  layer 
analysis  itself. 

The  streamline  routine  that  generates  data  for  use  in  the  boundary  lay¬ 
er  code  has  been  modified  to  correct  the  convergence  term.  The  convergence 
term  measures  how  quickly  the  surface  velocity  vectors  appear  to  converge. 
For  external  flow  over  an  axisymmetric  body  it  is  equal  to 

(l/r)(dr/dx) 


where  r  is  the  radius  of  the  body.  For  any  body  the  convergence  is 

^-1/vj  ^a(v  .  t)/atj 

where  t  is  the  vector  transverse  to  the  streamline  direction.  This  is  iden¬ 
tical  to  the  equation  for  a  body  of  revolution. 
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The  streamline  routine  at  one  time  made  the  assumption  that  the  surface 
potential  satisfied  Laplace's  equation  in  two  dimensions.  That  is, 

^av/asj  +  ^a(v  •  t)/atj  -  o 

Under  this  assumption  the  convergence  term  was  calculated  as 


[+i/vj  ^av/asj 

which  is  a  simple  calculation. 

Unfortunately,  no  physical  basis  exists  for  the  abova  assumption. 
The  simplest  illustration  of  this  is  the  flow  at  the  centorline  of  a  high 
aspect  ratio  wing.  The  surface  streamlines  have  zero  convergence  because 
they  remain  parallel;  however,  the  above  equation  will  produce  sizable 
values.  The  correction  of  this  error  is  even  more  important  to  internal 
flows  where  a  converging  duct  is  associated  with  Increasing  velocities. 

The  correct  general  formula  has  been  coded  into  the  streamline  routine 
and  the  computed  values  are  written  to  the  plot  file  for  display  as  the 
variable,  Rf,  under  the  streamline/two-dlmensional  option. 

Boundary  Laver  Methods 

The  Integral  method  is  an  adaptation  by  Curie  (17)  of  a  method  develop¬ 
ed  by  Thwaites  (18).  In  Thwaites'  method,  the  momentum  integral  equation, 

d^/dx  -  C^/2  -  (H  +  2)(fi/U)(dV/dx)  -  (tf/r)(dr/dx)  (3.1) 

is  written  in  the  form, 

(d/dx)(K/U)  -  L/U  -  (K/r)(dr/dx)/(dU/dx)  (3.2) 


where 


K  -  SV‘'(dU/dx) 

L  -  (i  -  K(H  +  2)) 

i  -  (^/U)(3U/ay)y_Q  (3.3) 

9  -  Momentum  Thickness 
H  -  Shape  Factor 
-  Skin  Friction 
u  -  Local  Velocity 

(l/r)(dr/dx)  -  a  measure  of  the  streamline  divergence  or  convergence. 
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In  Thwaites'  method,  the  divergence  term  is  not  considered,  although 
the  method  was  later  extended  to  Include  this  term  by  Rott  and  Crabtree 
(19). 


Thwaites  used  exact  solutions  to  a  variety  of  laminar  flows  to  deter¬ 
mine  the  relationship  between  L  and  K, 

L  -  0.^5  -  6K  (3.4) 

Curie  has  pointed  out  that  Eq.  (3.4)  is  not  adequate  in  flows  approaching 
separation,  and  he  has  suggested  an  extension  or  correction  giving 

L  -  0.45  -  6K  +  g(K,M)  (3.5) 

The  parameter,  /i,  is  a  function  of  both  the  pressure  gradient  and  the  curva¬ 
ture  or  second  derivative  of  velocity, 

/i  -  K2U(d2U/dx2)/(dU/dx)2  (3.6) 

Curie  rewrote  Eq.  (3.5)  in  the  form 

L  -  F  (K)  -  /iC  (K)  (3.7) 

o  o 

where  and  G  are  universal  functions  determined  from  a  series  of  exact 
solutions  to  laminar  flows  in  the  same  ws-  is  they  are  in  Eq.  (3.4).  After 
substitution  of  Eq.  (3.5)  into  Eq.  (3  2;,  and  with  subsequent  Integration, 
the  result  can  be  rearranged  Ir  the  form 

e  -  0.45i//(USr2)  f  r^d  2.22g)U5dx  (3.8) 

^  o 


This  equation  is  conveniently  solved  by  iteration,  g  initially  equal  to 
zero.  With  values  of  K  and  n  determined  in  the  first  iteration,  a  second 
iteration  is  carried  out  using  Eq.  (3.7).  At  each  step  in  the  calculation, 
the  local  skin  friction  coefficient,  C^,  and  the  shape  factor,  H,  can  be 
calculated  using  Eq.  (3.3).  The  local  skin  friction  coefficient  has  been 
defined  as 

Cf  -  (/i/p«U)i  (3.9) 

where  i  in  Eq.  (3.3)  is  determined  in  a  similar  manner  to  L  from  a  series  of 
known  solutions  to  give 

-  F^(K)k  -  G^(K)  (3.10) 

The  functions,  Fq,  F,  ,  Gq  and  Gj^  are  tabulated  in  the  computer  program. 
Calculations  begin  at  the  stagnation  point,  with  the  initial  momentum  thick¬ 
ness,  6,  given  as  a  function  of  K.  For  bluff  bodies,  K  takes  an  initial 


19 


o 


(3.11) 


r.0604l  ' 

[dU/dxJ 

The  calculation  proceeds  either  to  laminar  separation  or  to  the  end  of  the 
airfoil,  whichever  occurs  first.  The  calculated  boundary  layer  development 
is  then  interrogated  to  determine  if  transition,  laminar  separation  or 
forced  transition  (boundary  layer  tripping)  has  taken  place.  If  any  of 
these  phenomena  have  occurred,  the  downstream  flow  is  assumed  to  be  turbu¬ 
lent  . 


Boundary  Laver  Transition  and  Laminar  Separation 

Boundary  layer  transition  is  a  very  complex  phenomenon.  Reynolds  num¬ 
ber  is  a  controlling  parameter,  but  it  has  been  shown  that  the  Reynolds  num¬ 
ber  at  transition  can  be  increased  a  considerable  amount  by  careful  elimina¬ 
tion  of  disturbances.  At  very  low  Reynolds  numbers,  laminar  boundary  layers 
are  stable  to  small  disturbances.  However,  at  higher  Reynolds  numbers,  the 
boundary  layer  is  unstable,  and  small  disturbances  can  be  amplified.  Ampli¬ 
fication  of  these  disturbances  causes  the  flow  to  become  turbulent.  The 
point  at  which  flow  breakdown  occurs  depends  on  the  strength  and  dominant 
frequency  of  the  initial  disturbance.  Disturbances  may  be  due  to  free- 
stream  turbulence,  surface  roughness,  noise  or  vibration  of  the  surface.  As 
there  is  no  detailed  analysis  of  the  transition  process,  transition  predic¬ 
tion  is  accomplished  by  means  of  empirical  correlations.  Granville  (Ref. 
20)  has  developed  a  procedure  based  on  the  determination  of  the  neutral  sta¬ 
bility  point  and  the  transition  point.  The  neutral  stability  point  is  de¬ 
fined  as  that  point  downstream  from  which  small  disturbances  are  amplified 
within  the  boundary  layer.  It  is  this  amplification  of  small  disturbances 
chat  ultimately  leads  Co  transition.  The  neutral  stability  point  is  reached 
when  the  Reynolds  number,  based  on  the  local  momentum  thickness  and  local 
flow  properties,  attains  some  critical  value,  R^j^  Schlichting  and 
Ulrich  (21)  have  shown  that  R. .  can  be  correlated^ w?fh  the  local  pres¬ 
sure  gradient  parameter,  K  -  (dU/ds)  .  Correlations  by  Smith  (22)  and 
others  have  been  reduced  to  analytical  form  as  follows. 


Instability  Curves 

K  -  -0.4709  +  0.11066  InR,  -  0.0058591  In^R 


6 


and 


for  0  <  R.,  <  650 

flins  - 


K  -  0.69412  -  0.23992  lnR„  +  0.0205  In^R 


B 


for  650  <  R 


ins 


(3.12a) 


(3.12b) 
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If  for  a  given  R,  the  pressure  gradient  parameter,  K,  as  calculated  by  Eq. 
(3.12),  is  greater  than  determined  by  the  boundary  layer  development,  the 
flow  has  passed  from  a  stable  to  an  unstable  region.  Once  the  flow  passes 
into  the  unstable  region,  the  transition  process  begins;  and  Granville  has 
been  able  to  show  that  a  correlation  similar  to  the  instability  process  can 
be  used  to  determine  the  transition  point.  He  formed  an  average  pressure 

gradient  parameter,  K,  defined  as 


K 


r^trans 

J^ins 


K  ds 


trans 


ms 


(3.13) 


which  correlated  reasonably  well  with  with  the  momentum  thickness  Reynolds 
number  at  transition,  R.  This  correlation  is  presented  in  analytical 
form  as  follows.  » trans 


Transition  Curves 


K  -  -  0.0925  +  7.0  X  10  R 


B 


K  -  -  0.12571  +  1.14286  x  lO'^  R 


8 


for  750  <  R.^  <  1,100 

5 trans  - 


and 


(3.14a) 


(3.14b) 


K  -  1.59381  -  0.45543  InR  +  0.032534  In^R. 

r  V 


(3.14c) 


for  1,100  <  R 


$ trans 


When  the  K  calculated  by  Eq.  (3.14)  for  a  given  R^  is  greater  than  the  value 
determined  from  the  boundary  layer  development,  transition  is  predicted. 

With  transition  predicted,  initial  values  of  the  momentum  thickness,  8, 
and  shape  factor,  H,  are  required  to  start  the  turbulent  boundary  layer  cal¬ 
culation.  Because  the  boundary  layer  growth  is  continuous,  the  momentum 
thickness  at  transition  is  used  as  the  initial  turbulent  momentvun  thickness. 
Since  the  shape  factor  varies  from  values  greater  than  2-0  to  less  than  1.5 
through  the  transition  region,  an  empirical  expression  is  used  to  determine 
the  initial  turbulent  shape  factor.  The  empirical  relation  between  H  and 

R.  was  determined  from  the  data  obtained  by  Coles  (23) : 

S trans 
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t  Log 


]  .4754 


+  0.9698 


10  ‘'tftrans 


(3.15) 


In  many  cases,  the  pressure  gradient  is  of  sufficient  strength  to  sepa¬ 
rate  the  laminar  boundary  layer  prior  to  transition.  Except  in  extreme 
cases,  the  boundary  layer  will  then  reattach,  usually  as  a  turbulent  bound¬ 
ary  layer.  Researchers  been  able  to  analyze  this  phenomenon  (24),  but  the 
procedure  is  extremely  complicated  and  cumbersome;  consequently,  empirical 
relationships  are  required.  From  the  measurements  of  Caster  (25)  and  oth¬ 
ers,  a  correlation  is  formed  which  is  capable  of  predicting  both  the  occur¬ 
rence  of  a  separation  and  later  the  reattachment  as  a  turbulent  boundary 
layer  or  the  catastrophic  separation.  The  correlation  is  of  the  form, 

K  -  0.0227  -  0.0007575  -  0.000001157  R^2  (3.16a) 
for  Rg  >  125 


and 


K  -  - .09 


for  R^  <  125. 


(3.16b) 


The  value  of  K  determined  by  the  boundary  layer  development  must  be 
less  than  that  calculated  by  Eq.  (3.16)  before  separation  without 
reattachment  is  predicted.  If  reattachment  is  predicted,  the  turbulent 
boundary  layer  calculation  is  initiated  using  the  momentum  thickness 
calculated  at  the  separation  point. 


Turbulent  Method 


Methods  for  the  calculation  of  turbulent  boundary  layers  in  two  dimen¬ 
sions  have  been  developed  by  many  investigators.  These  methods  were  re¬ 
viewed  at  a  conference  held  in  1968  at  Stanford  University  (26).  One  of  the 
methods,  an  integral  method  by  Nash  and  Hicks  (27),  compared  very  favorably 
with  the  more  complex  finite  difference  methods.  The  method  remains  (White 
(28))  an  excellent  approach  for  application  to  the  current  problem,  both  in 
terms  of  accuracy  and  speed. 

The  Nash-Hicks  method  is  based  on  momentum  and  moment  of  momentum  equa¬ 
tions  coupled  with  a  skin  friction  law  derived  from  Coles'  velocity  profile 
family  (29).  An  additional  equation  is  obtained  by  relating  the  shear 
stress  integral  to  its  equilibrium  value  using  a  simple  first-order  differ¬ 
ential  equation.  The  equations  have  been  derived  in  (30)  and  (31),  and  are 
repeated  herein  for  completeness. 
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A  family  of  integral  equations  can  be  derived  taking  higher  moments  of 
the  equation  of  motion.  The  resulting  equations  can  be  expressed  as 


1  Sm 

r  dy 


Slail 

dx 


a.  Ul"'  dU  1  r*  a  Sr . 

y  “y  -  fn  ;  J  y  ^  <*y  <’■*’> 


where 

a  -  0  gives  the  momentum  integral  equation,  and 
a  -  1  gives  the  moment  of  momentum  integral  equation 


The  velocity  distribution  across  the  boundary  layer  can  be  represented  by 
Coles'  velocity  profile  family,  given  by 


u 


u  u 

^  {ln(y  ^)  +  C)  + 
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{1 


cos 


S  ’ 


(3.18) 


where 


K 


c 


1  /2 

friction  velocity,  (r ^p)  '  , 

free  parameter  having  units  of  velocity, 

boundary  layer  thickness , 

.41,  and 
2.05. 


Substitution  of  this  equation  into  Eq.  (3.17)  results  in  two  equations  of 
the  form; 


+  B  f  -  D 

dx  (  dx  dx  r  dx 


(3.19) 


A  third  equation  of  the  same  form  is  obtained  by  evaluating  Eq.  (3.18)  at  y 
-  S  followed  by  the  differential  with  respect  to  x.  The  parameter,  ♦,  is 
represented  by  the  following  relations: 
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(i)  a  -  0 


-1  ^ 

♦  -  f  u 

6  T 


(li)  a  -  1 


f  " 


(iii)  a  —  « 


*  -  0 


(3.20) 


rS  ^ 

The  shear  stress  integral,  I  p  ^y>  appearing  in  Eq.  (3.20)  was  evaluated  by 
Nash  and  Hicks  using  an  equation  of  the  form: 


where 


(3.21) 


-  -  2 

1/2  pU  S  Jo 

The  equilibrium  value  of  ,  (C^g^) .  was  determined  by  Nash  and  Macdonald 
(30)  from  measured  shear  stress  distributions  giving 


C  -  .025  (1  -  i)* 
req  H 


(3.22) 


where  H  is  the  local  shape  factor.  Equation  (3.22)  can  be  expressed  in 
terms  of  the  parameters,  u  ,  u„  and  U,  by  evaluation  of  the  integral  rela¬ 
tions  used  to  define  H;  that  is, 


-  - 


(1  -  P  dy 


I S 


(3.23) 


(I  -  dy 


with  the  aid  of  Eq.  (3.18).  The  equations  for  the  four  unknowns,  u  ,  u„,  S, 
and  ,  can  be  integrated  once  the  pressure  distribution,  U(x) ,  and^  the 
streamline  divergence,  l/r(dr/dx),  are  prescribed.  Initial  conditions  are 
obtained  from  the  transition  analysis.  In  this  case,  initial  values  of  the 
momentum  thickness,  $,  and  the  shape  factor,  H,  are  known.  Initial  values 
of  u^  ,  u^  and  6  can  be  determined  from  the  known  8  and  H  using  Eq.  (3.18). 
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The  starting  value  for  C  is  obtained  by  making  the  assumption  that  in  the 

region  of  transition,  C  =  C 

r  req 

The  accuracy  of  the  calculation  method  is  demonstrated  in  Figure  3.1. 
This  figure  shows  a  comparison  between  measured  and  calculated  boundary  lay¬ 
er  developments  along  a  streamline  on  the  U.S.  Airship  Akron.  A  further 
comparison  is  shown  in  Figure  3.2  for  the  case  of  a  boundary  layer  in  a 
strong,  adverse  pressure  gradient  approaching  separation.  Particular  refer¬ 
ence  should  be  made  to  the  good  agreement  between  calculated  and  measured 
skin  friction  coefficients. 


3 . 2  Applications 

Several  test  cases  have  been  constructed  to  verify  the  working  of  the 
code  under  conditions  of  engineering  interest.  Cases  involving  radially 
nonuniform  flow  through  propeller,  hub  and  duct  assemblies  have  been  run 
successfully.  The  propeller  chosen  had  a  skew  of  72®  with  a  very  small 
clearance  of  O.OIR  between  the  duct  wall  and  the  propeller  tip.  Multiple 
stator  blades  were  constructed  using  the  NACA  0012  profile.  These  were  po¬ 
sitioned  forward  of  the  propeller.  Separate  cases  were  run  with  and  without 
the  stators.  The  flow  with  stator,  hub  and  duct  was  calculated  with  a  ve¬ 
locity  survey  at  the  grid  points  of  a  three-dimensional  mesh.  The  calcula¬ 
ted  velocity  was  used  in  a  separate  run  as  a  nonuniform  inflow  around  the 
propeller,  hub  and  duct  assemblies.  Provision  was  made  for  the  user  to  pro¬ 
vide  the  nonuniform  inflow  data  either  by  means  of  a  separate  tape  file  or 
as  part  of  the  input  file.  The  latter  option  is  easier  when  radially  non- 
uniform  inflow  is  present  or  when  the  volume  of  inflow  data  is  not  very 
large.  The  former  mode  of  specification  requires  some  user  manipulation  of 
input  such  as  providing  the  data  in  the  auxiliary  tape  unit  in  the  desired 
format.  For  minimum  user  manipulation,  a  direct  solution  with  all  the  com¬ 
ponent  parts  in  place  would,  of  course,  be  preferable. 

The  number  of  panels  used  on  the  propeller  was  320,  and  was  not  regard¬ 
ed  as  sufficient  for  accurate  Cp  calculation.  Five  blades  were  used  in  all 
propeller  calculations.  Hub  and  duct  were  input  as  axisymmetric  components 
over  a  basic  region  of  angular  width,  360/5  -  72®.  When  a  different  number 
of  blades  was  used,  an  appropriate  change  must  be  made  so  that  with  NBLAlJE 
rotational  reflections,  a  full  body  was  simulated.  The  cases  documented  are 
as  follows. 

1.  Five-bladed  propeller  in  axisymmetric  duct  and  hub.  The  hub  was  re¬ 
garded  as  fixed.  It  may  be  changed  to  rotating  by  a  single  input 
change.  The  entire  assembly  was  subject  to  nonuniform  radial  flow  as 
specified  above. 

2.  Five-bladed  stator  (NACA  0012  section)  in  axisymmetric  duct  and  hub 
subject  to  radially  nonuniform  flow. 
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PREDICTED 


Figures  3.3(a)  and  (b)  are  two  views  of  the  propeller,  hub  and  duct  as¬ 
sembly.  As  mentioned  earlier,  the  duct  and  hub  are  defined  over  a  basic  re¬ 
gion  which  is  reflected  NBLADE  times  to  provide  complete  definition  over 
360*.  Figures  3.4(a),  (b)  and  (c)  show  the  sectional  Cp  distribution  at  the 
three  radial  locations.  The  Cp  distribution  without  the  duct  and  the  hub 
(Figure  3.5)  is  also  shown  to  highlight  the  difference.  The  solution  near 
the  hub  is  more  strongly  influenced  compared  to  the  tip  region.  Figures 
3.6(a)  and  (b)  are  sectional  Cp  plots  along  a  longitudinal  (constant  Y)  sta¬ 
tion  and  a  lateral  station  (constant  X).  Figure  3.7(a)  and  (b)  show  the  im- 
final  calculation  without  a  stator. 

Figures  3.8(a)  and  (b)  are  views  of  the  five-bladed  stator,  duct  and 
hub.  Figure  3.9  shows  the  sectional  Cp  plot  on  the  primary  stator.  The  Cp 
is  symmetric  as  the  stators  are  carrying  zero  load.  The  inflow  is  radially 
nonuniform  as  before,  but  preserves  symmetry.  Figure  3.10  shows  contour 
plots  of  a  velocity  scan.  Care  is  necessary  in  interpreting  the  velocity 
scan  from  a  rotating  propeller.  The  flow  is  time -dependent  and  the  value  of 
velocity  at  a  point  is  what  an  observer  rotating  with  the  point  at  the  an¬ 
gular  speed  of  the  propeller  will  see  as  a  steady  velocity.  An  azimuthal 
averaging  scheme  has  been  added  which  provides  a  time -averaged  velocity  cal¬ 
culation  in  addition  to  the  time  history  of  the  velocity  components.  The 
velocity  scan  models  in  an  approximate  but  realistic  way  the  influence  of 
the  stator,  hub  and  duct  assembly  on  the  surrounding  flows. 
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(b)  Oblique  View  Showing  Shroud  Paneling 
Figure  3.3.  Concluded. 


NO  HUB.  NO  DUCT 


essure  Distribution  on  Propeller  with  Shroud  and  Hub. 


Figure  3.5.  Pressure  Distribution  in  Radially  Non-Uniform  Flow. 


NO  STATOR 


Propeller,  Hub  and  Shroud  Assembly  In  the  Non-Uniform  Flow  Induced  by  Stator,  Hub  and  Shroud. 


NO  STATOR 


Figure  3.7  .  Concluded. 
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Figure  3.8.  Concluded. 


OFF  BODV  DATA  BLOCK  1  I -PLANE 


Figure  3.10.  Velocity  Scan  Data  Showing  Flow  Non-Uniformity. 


4.0  ISOLATED  HULL  ANALYSIS 


An  accurate  prediction  of  the  flow  at  the  stern  is  of  great  practical 
interest  in  the  determination  of  propeller  hull  interaction.  The  improve¬ 
ments  to  the  boundary  layer  and  the  streamline  analysis  reported  in  the  pre¬ 
vious  chapter  aid  the  correct  modeling  of  hull  flow.  There  is,  however,  an 
additional  problem  associated  with  the  analysis  of  the  viscous  flow  over  the 
stern  arising  from  the  failure  of  conventional  boundary  layer  calculation 
methods  to  account  for  thick  boundary  layers.  In  the  literature  this  prob¬ 
lem  has  been  addressed  in  two  general  directions.  One  of  these  involves 
generalization  and  extension  of  the  thin  boundary  layer  methods,  of  both 
integral  and  differential  types,  to  include  such  factors  as  changes  in  coor¬ 
dinate  metrics  normal  to  the  surface,  normal  variation  of  pressure,  and  in¬ 
teraction  between  viscous  and  inviscid  flow  regions.  The  alternative  ap¬ 
proach  to  the  calculation  of  stern  flows  involves  numerical  solutions  to  the 
complete  Reynolds  -  averaged  Navier-Stokes  equations. 

In  the  present  work  both  of  these  approaches  have  been  applied  in  a 
limited  way.  MPROP  has  been  applied  to  -alculate  the  potential  flow  about 
the  hull  form.  The  approximate  bouu<!ary  layer  profiles  calculated  by  MPROP 
have  been  used  in  an  iterative  scheme  to  model  interaction  between  viscid 
and  inviscid  flows.  The  Reynolds -averaged  Navier-Stokes  solver  developed  by 
Patel  et  al.  (32)  was  employed  to  calculate  the  viscous  velocity  profile  in 
the  plane  of  the  propeller.  Both  approaches  are  briefly  described  below. 


1  Potential  Flow 

The  boundary  layer  characteristics  of  the  Akron  Airship  (Figures 
4.1(a),  (b)  and  4.2(a),  (b))  were  calculated  before  and  after  the  changes  to 
MPROP.  The  accuracy  of  the  calculation  method  is  demonstrated  in  Figure 
4.3.  This  figure  shows  a  comparison  between  measured  and  calculated  bound¬ 
ary  layer  developments  along  a  streamline  on  the  U.S.  Airship  Akron. 

The  new  and  old  calculations  are  compared  at  zero  angle  of  attack  and  a 
Reynolds  number  of  4.56  million.  As  seen  in  Fig.  4.1,  the  momentum  thick¬ 
ness  is  now  considerably  larger  near  the  tail.  The  prediction  was  25%  below 
experiment;  it  is  now  10%  above  the  last  experimental  datum.  The  boundary 
layer  thickness  was  also  25%  below  experiment  with  the  wrong  convergence 
term.  Fig.  4.2(a).  It  is  now  10%  below.  The  skin  friction  and  shape  factor 
were  only  minutely  affected  by  the  change.  Two  comparisons  with  experiment 
are  shown  in  Fig.  4.3;  one  near  the  midpoint  of  the  hull,  the  other  near 
the  tail.  It  is  worth  noting  that  the  good  agreement  would  not  be  possible 
unless  both  the  external  velocity  and  boundary  layer  thickness  matched  that 
of  experiment. 
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Inside  Boundary  Layer  Correlated  with  Experiment. 


4.2  Viscous  Flow 


Axlsynunetric  bodies  were  chosen  for  this  investigation  in  order  to  di¬ 
rectly  focus  on  the  interaction  between  a  propeller  and  a  thick  stern  bound¬ 
ary  layer.  Experimental  results  for  a  set  of  body  of  revolution  models  are 
available  (33).  Figure  4.4  shows  the  geometry  of  the  axisymmetric  after¬ 
bodies.  For  the  present  investigation.  Analytical  Methods  has  considered 
two  different  codes  to  calculate  the  flow  around  the  hull.  These  are  the 
Reynolds -averaged  Navier-Stokes  solver  (RANSTERN)  and  ARC3D.  RANSTERN  is  a 
numerical  method  developed  by  Chen  and  Patel  (34)  for  the  solution  of  the 
fully-elliptic  Reynolds -averaged  Navier-Stokes  equations.  ARC3D,  which  was 
originally  developed  by  Steger  at  NASA  Ames  Research  Center  (35)  is  based  on 
the  implicit  approximate  factorization  algorithm  of  Beam  and  Uarming  (36) 
for  unsteady  compressible  flow. 

Figures  4.5(a)  and  (b)  show  the  grid  employed  to  calculate  the  flow  on 
Afterbody  1.  The  pressure  distribution  on  Afterbody  1  at  a  Reynolds  number 
of  6.6  million  is  compared  with  experiment.  Both  ARC3D  and  RANSTERN  cor¬ 
rectly  predict  the  pressure  on  the  body  surface.  The  second  pressure  peak 
is  more  accurately  predicted  by  RANSTERN.  Further  analysis  for  the  calcula¬ 
tion  of  wake  flow  was  made  using  RANSTERN. 
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Axlsymmetrlc  Afterbodies 


CTX 


Pressure  Distribution  on  Afterbody  1  at  Re  ■  6.6  Million. 


5.0  PROPELLER -HULL  INTERACTION 


The  complete  problem  of  propeller-hull  interaction  without  restrictive 
assumptions  on  geometry,  the  nature  of  viscous  effects  and  types  of  cavita¬ 
tion  is  both  difficult  and  subject  to  prohibitive  computational  costs.  Many 
of  the  physical  processes  involved  in  propeller-hull  interaction  are  locally 
dominated  by  viscosity.  However,  it  is  still  appropriate  to  represent  the 
flow  field  due  to  the  propeller  action  by  means  of  a  velocity  potential  sa¬ 
tisfying  the  Laplace  equation.  Of  the  different  computational  approaches 
available  to  solve  the  resulting  equations,  the  Green's  function  surface 
singularity  method  coupled  with  special  routines  for  nonlinear  effects  has 
the  promise  of  success. 


5 . 1  Technical  Approach 

Consider  a  propeller  operating  in  the  wake  of  a  long  hull.  The  radius 
of  the  hull  wake  is  assumed  to  be  comparable  in  size  to  the  propeller  radi¬ 
us.  In  other  words,  the  propeller  is  completely  immersed  in  the  boundary 
layer  of  the  hull.  For  the  purposes  of  the  analysis  the  following  two  vis¬ 
cous  assumptions  are  made. 

(1)  The  viscc..s  hull  flow  at  the  propeller  plane  is  calculated  by  a 
three-dimensional  finite  difference  code  by  neglecting  the  presence 
of  the  propeller. 

(2)  The  unsteady  propeller  flow  is  modeled  as  a  potential  flow  problem 
with  nonuniform  inflow  describing  the  hull  wake. 

By  the  above  two  assumptions  the  overall  problem  of  propeller-hull  in¬ 
teraction  becomes  decoupled  into  two  constituent  parts.  The  first  part  is 
entirely  viscous  and  three  dimensional.  Finite-difference  codes  solving 
Navier- Stokes  equations  are  appropriate  for  the  viscous  analysis.  The  vis¬ 
cous  analysis  provides  the  wake  fractions  in  the  plane  of  the  propeller. 
The  propeller  field  is  constructed  by  the  method  of  surface  singularity 
distribution  satisfying  an  implicit  boundary  condition  of  the  exit  plane  of 
the  viscous  domain.  Inviscid  flow  techniques  are  applicable  because  the 
vorticity  contained  in  the  hull  wake  is  neglected. 

The  scope  of  the  present  investigation  is  limited  to  calculating  the 
flow  field  around  a  propeller  operating  in  the  presence  of  a  thick  and  pos¬ 
sibly  separated  turbulent  hull  wake.  The  focus  here  is  on  the  unsteady 
pressure  distribution  on  the  propeller  surface  and  is  caused  by  the  azi¬ 
muthal  variation  of  Inflow  in  the  plane  of  the  propeller.  This  problem  has 
been  modeled  in  two  different  ways.  Firstly,  MPROP  has  been  modified  to 
calculate  a  series  of  solutions  as  a  perturbation  over  a  basic  uniform  flow 
solution.  Secondly,  the  entire  problem  has  been  analyzed  by  USAERO,  an  un¬ 
steady  code  marched  in  time  through  several  revolutions  of  the  propeller. 
In  the  following  paragraphs,  a  brief  description  of  the  two  methods  are 
given. 
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5.2 


Frequency  Band  Analysis 

Consider  an  assemblage  of  several  three-dimensional  closed  bodies  with 
an  outer  boundary,  S,  moving  through  a  large  mass  of  constant  density  fluid. 
Let  W  denote  the  wake  surface  of  bodies  with  non-zero  circulation.  The  ve¬ 
locity  potential  across  the  wake  is  discontinuous.  In  the  region  excluding 
the  potential  barrier  of  the  wake,  the  flow  is  non- circulatory  and  the  ve¬ 
locity  potential  at  any  point,  P,  in  the  flow  (following  Lamb's  (15))  is: 


^  lx  ,y  ,z  ,tl  - 

pI  P  P  J 


1_ 

4ir 


1  M 

r  3n 


dS  + 
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II 

s+w 
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(5.1) 


where 


1 

n - p — - p — - p 

y  X  -  X,  +  y  -  Yi  +  z  -  z. 

Ip  iJ  I  p  iJ  Ip  iJ 


(5.2) 


and  X.  ,yj^  and  represent  a  point  on  the  body  surface  corresponding  to  the 
variaole^of  integration.  The  integration  must  extend  over  all  points  on  the 
body  and  wake  surface.  4  is  the  disturbance  potential  due  to  body  motion 
and  vanishes  at  infinity.  When  point  P  is  on  the  surface,  the  locally  sin¬ 
gular  kernel,  1/r,  must  be  replaced  by  the  appropriate  limiting  value  and 
should  be  excluded  from  the  region  of  integration. 

Discretization  of  S  into  N  panels,  assuming  that  ^  is  constant  on  each 
of  the  panels  and  d^/dn  takes  the  value  of  prescribed  normal  velocity,  pro¬ 
vides  the  following  matrix  equation  for  the  unknown  velocity  potential  on 

the  surface,  ^  . 

s 


[A]  [^3]  -  [^j,] 

Elements  of  matrix  A  are  geometric  quantities  defined  by 

J|  a/an  (1/r)  dS 

and  is  the  velocity  potential  induced  at  a  control  point  by  the  sum  of 

all  constant  strength  sources  of  known  strength.  The  former  may  be  recog¬ 
nized  as  the  velocity  potential  induced  at  a  control  point  by  a  unit 
strength  doublet  whose  axis  is  normal  to  the  panel  surface.  The  above  ma¬ 
trix  equation  is  the  basic  equation. 

The  inclusion  of  a  spatially  non-uniform  velocity  distribution  result¬ 
ing  from  hull  wake  may  now  be  considered.  Th;>  flow  in  the  steady  state  must 
repeat  itself  every  2»r/N  degrees  of  rotation  even  when  spatially  non-uniform 
inflow  from  the  hull  wake  is  present.  By  comparison  with  the  steady  propel¬ 
ler  problem,  the  major  change  is  circumferential  variation  of  velocity  nor¬ 
mal  to  the  propeller.  As  the  blade  rotates  it  encounters  a  change  in  the 
normal  velocity  as  a  function  of  $,  the  angular  position  of  the  blade.  The 
general  problem  would  appear  to  rule  out  the  possibility  of  employing  a 
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method  based  on  steady  flow  analysis.  However,  there  appears  to  be  a  way 
whereby  the  steady  flow  procedure  can  be  applied  to  model  general  non-uni¬ 
form  flow.  The  application  of  that  procedure  will  be  outlined  below.  It  is 
assumed  that  non-uniform  flow  is  prescribed. 


Let  the  non-uniform  inflow  be  prescribed  as  a  function  of  r  and  $; 
l.e. .  V  -  f(r.tf). 


2n  R 

Let  the  propeller  rotate  at  n  revolutions  per  second.  J  J  V  dr  dfl  - 

Q,  say.  Now,  define  V  -  Q/xR^ ,  where  R  is  the  propeller  rad?us .  Now  the 
prescribed  non-uniform  inflow  can  be  represented  as 


V  -  Vq  -  cos  5  +  82  cos  26  +  a^  cos  36  +  ... 

+  b^  sin  +  b2  sin  26  +  b^  sin  36  (5.3) 


Typically,  a^^,  a^,  ....  b^^,  b2,  etc.  are  functions  of  r  and  small  in  magni¬ 
tude  compared^ to  v^. 

A  solution  is  next  sought  for  the  propeller  in  non-uniform  inflow  as 
follows . 


^  ^  ^  cos  cos  26  +  cos  36  ... 

q  Ic  2c  3c 

+  sin  6  +  sin  26  +  sin  36  (5.4) 

Is  2s  3s  ' 

4  is  the  unknown  doublet  potential  on  a  panel.  ^  is  the  potential  value 
based  on  steady  propeller  flow  at  an  advance  coeffic?ent,  J  -  V_/nD.  ^ic’ 
. ^2^  .  ^2s'  unknowns  to  be  determined.  Hire  tne  relation¬ 
ship  between  Blade-fixed  coordinate  system  and  space-fixed  coordinate  system 
may  be  noted.  The  onset  flow  is  defined  by  Eq.  (5.3)  above  in  which  6  is 
measured  counterclockwise  from  the  z-axis.  A  panel  on  a  propeller  rotating 
at  w  radians  per  second  will  negotiate  an  onset  velocity  corresponding  to  ^ 
-  wt,  where  ^  is  the  angle  made  by  the  panel  control  point  relative  to  the 
blade  fixed  axis;  i.e.,  6  -  ^  -  wt. 

The  coefficients  of  and  are  to  be  determined  from  a  knowledge 
of  source  strengths  required  to  account  for  the  circumferential  variation  of 
onset  flow.  For  a  simple  example,  assume  the  onset  flow  is  given  by 

V  -  Vq  +  cos  6 
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A  solution  is  sought  in  the  form 


Now,  on  the  blade  surface, 


^  cos 


3^/3n  -  V  •  n 


or 


r  1 

_ +  _  cos  -  V  +  a. 

^  fln  ;)t,  '  W  i- 


cos  0 


(5.5) 


dn  dn 


Thus,  d^,/dn  -  a,  will  satisfy  the  boundary  conditions  on  normal  velocity. 
Based  on  Eq.  (5.1)7 

♦ic  -  •  b  J1  i  *1  <**  *  h  IJ  *ic  t  <i) 


s+w 


The  above  equation  is  exactly  equivalent  to  a  doublet  solution  in  response 
to  a  panel  source  strength  of  a...  Though  a.  is  a  constant,  it  actually 
varies  as  a  function  of  panel  angular  position,  0  -  wt.  In  other  words, 
with  the  source  strength,  a.  being  knovm,  it  is  possible  to  solve  for  in 
a  manner  similar  to  that  of ^steady  flow  analysis.  ^ 

An  important  aspect  of  unsteady  flow  is  the  shedding  of  variable 
strength  vorticity.  Ideally,  a  harmonic  wake  (i.e.,  a  wake  in  which  the 
strength  varies  as  cos  (x/i)  where  x  is  the  distance  along  the  wake  and  i  is 
the  net  distance  travelled  by  the  wake  in  one  revolution  of  the  blade)  must 
be  modeled.  In  practise,  these  integrals  can  be  constructed  in  a  closed 
form.  At  this  point  a  simplifying  assumption  is  introduced- -the  response  of 
the  propeller  perturbation  potential  to  the  circumferential  velocity  com¬ 
ponent  is  assumed  to  be  in  phase .  The  phase  difference  between  the  onset 
flow  and  wake  strength  has  been  neglected.  Under  this  assumption  the  prob¬ 
lem  becomes  similar  to  that  of  finding  the  solution  for  a  steady  onset  flow 
of  magnitude  a^.  Now,  the  unsteady  pressure,  d^/dt,  can  be  evaluated,  re¬ 
membering  that  0  -  wt  -  Ip. 

dp/dt  -  -aw  cos  (wt  -  0^^) 

for  the  simple  example. 
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5 . 3  Time -Dependent  Analysis 
5.3.1  Background 


The  time -dependent  analysis  of  a  propeller  in  nonuniform  flow  was  com¬ 
puted  using  an  existing  program,  USAERO.  This  code  was  developed  at  Analy¬ 
tical  Methods,  Inc.  originally  for  helicopter  blade  tip  analyses  (37).  Lat¬ 
er,  the  method  was  extended  to  treat  more  general  problems.  The  current 
version  calculates  the  unsteady  subsonic  flow  characteristics  for  multiple 
bodies  in  general  motion.  The  motions  may  be  prescribed  or  they  can  be  cal¬ 
culated  internally  using  a  six-degree-of -freedom  flight  path  integration 
routine  which  is  an  optional  module. 

Based  on  a  time-stepping  approach,  USAERO  provides  a  very  general  capa¬ 
bility  for  applications  to  such  transient  problems  as  helicopter  rotor/fuse¬ 
lage  interactions  (38),  vehicle-on-vehicle  interactions,  propeller/stator 
problems,  maneuvering  vehicle  (39)  and  store  release  problems.  Output  from 
the  program  includes  time  histories  of  integrated  force  and  moment  charac¬ 
teristics  for  selected  parts  of  a  CONFIGORATim  together  with  the  corre¬ 
sponding  velocity  and  pressure  distributions.  These  results,  together  with 
the  three-dimensional  geometry  of  the  CONFIGURATION  (including  the  generated 
wakes) ,  may  be  displayed  graphically  step  by  step  using  the  interactive  pro¬ 
gram,  OMNIPLOT,  on  a  graphics  workstation  or  in  an  "animated"  sequence  using 
interactive  program,  0MNI3D,  on  a  Silicon  Graphics  IRIS  4-D  series  worksta¬ 
tion. 


5.3.2  Method  (Xitline 

The  basis  of  the  method  is  a  time-stepping,  surface  singularity  method 
which  uses  quadrilateral  PANEL2  of  uniformly  distributed  doublets  and 
sources.  The  surface  integrals  in  Green’s  theorem  are  evaluated  in  a  piece- 
wise  manner  over  each  PANEL  to  form  panel  influence  coefficients.  These  are 
evaluated  for  each  PANEL  acting  at  the  central  control  points  on  all  the 
surface  PANELS,  thus  forming  a  matrix  of  influence  coefficients.  Usually 
the  source  values  are  determined  at  the  start  of  each  time  step  according  to 
the  local  velocity  component  normal  to  the  PANEL  surface  (the  source  values 
include  terms  for  normal  inflow/outflow  as  well  as  boundary  layer  displace¬ 
ment)  .  The  doublet  values  are  then  solved  using  the  matrix  equations. 

There  is  also  provision  to  solve  the  source  values  (i.e.,  normal  veloc¬ 
ity)  given  the  panel  doublet  strengths  (potential) .  This  is  done  on  a  PATOl 
by  PATCH  basis  (a  PATCH  is  a  collection  of  panels  in  a  rectangular  array) . 
Thus  mixed  boundary  condition  problems  may  be  treated.  Some  PATCHES  may 
have  user  input  values  of  both  source  and  doublet  for  special  flow  modeling 
options . 

Extensive  component  copying  facilities  are  provided  at  the  input  level. 
These  not  only  provide  a  convenient  way  of  describing  the  multiple  blades  of 
a  propeller,  but  also  provide  an  option  for  blades  to  copy  earlier  doublet 
solutions  from  a  "primary"  blade  in  situations  where  the  solution  repeats  as 
blades  pass  the  same  azimuthal  station.  This  capability  is  applicable  in 
the  case  of  the  propeller  in  nonuniform  flow  provided  the  nonuniformity  is 
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steady  in  the  ground- fixed  frame.  This  reduces  the  number  of  unknowns,  and 
hence  computing  time,  considerably. 

The  tangential  velocity  at  the  solid  surfaces  of  the  CONFIGURATION  is 
obtained  from  the  doublet  gradient.  The  normal  component  (which  can  be  non¬ 
zero  for  inflow/outflow  and  transpiration  modeling)  is  provided  by  the 
source  value.  Off -body  velocity  perturbations  are  evaluated  by  summing  all 
doublet  and  source  singularity  contributions.  These  velocities  are  used  for 
general  flow  field  information  as  well  as  for  wake  point  convection.  The 
WAKES  grow  with  each  time  step,  with  new  wake  points  being  propagated  from 
wake -shedding  lines,  while  all  the  previous  wake  points  are  convected  at  the 
local  fluid  velocity. 


5.3.3  Method  Formulation 
Basic  Equations 

Consider  the  CONFIGURATION  moving  with  velocity,  Vg,  through  an  un¬ 
bounded  fluid  initially  at  rest.  Figure  5.1.  The  basic  assumptions  are  that 
the  effects  of  viscosity  are  largely  confined  to  thin  boundary  layers  on  the 
CONFIGURATION  surface  and  that  wake  vorticity  is  essentially  concentrated  in 
thin,  free -shear  layers  and  discrete  vortex  filaments.  The  majority  of  the 
flow  is,  therefore,  regarded  as  inviscid,  irrotational  and  incompressible. 
Laplace’s  equation  can  then  be  applied, 

7’^  -  0  (5.6) 

The  convention  adopted  here  is  that  the  perturbation  velocity  is  the  nega¬ 
tive  gradient  of 


V  -  -7^ 


(5.7) 


Green's  Theorem  is  applied  next.  Note  that  with  7^^  -  0  the  volume 
integral  disappears.  The  flow  is  therefore  uniquely  determined  by  surface 
integrals  of  ^  and  its  normal  derivative  over  the  surface  of  the  CONFIGURA¬ 
TION  and  its  WAKE.  Thus  the  velocity  potential,  ^  for  a  point,  P,  on  the 
wetted  side  of  the  surface  is  ^ 


I  I  *  n  .  7  [1)  dS  +  ^  I  I  1  n  .  dS 

S-P  s 

*  i;  n  ’  (t] 


(5.8) 


where  n  is  the  outward  normal  from  the  surface  and  r  is  the  length  of  the 
vector  from  the  surface  element,  dS,  to  the  point,  P.  S-P  signifies  that 
the  point,  P,  is  excluded  from  the  surface  Integral- -  the  limiting  process 
for  the  singular  point  when  r  -*  0  yields  the  local  contribution,  ^p/2. 
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Figure  5.1.  General  Arrangement. 


The  first  Integral  In  Eq.  (5.8)  is  the  contribution  from  a  surface  dis¬ 
tribution  of  normal  doublets  of  strength, 

/i  -  ^/4ir  (5.9) 


The  second  integral  is  the  contribution  from  a  surface  distribution  of 
sources  of  strength, 


a 


n  «  7<A 
Air 


(5.10) 


The  third  integral  in  Eq.  (5.8)  is  the  contribution  from  the  wake  sur¬ 
face.  Here,  the  upper  and  lower  surfaces  have  been  combined,  taking  the  up¬ 
ward  facing  normal.  The  upper  and  lower  potentials  have  been  combined 
yielding  a  wake  doublet  distribution  of  strength, 


-  ^L> 

^  "  Air 


(5.11) 


This  is  the  potential  Jump  across  the  wake.  In  combining  the  upper  and  low¬ 
er  surfaces  the  source  term  has  been  discarded.  Implying  that  there  is  no 
normal  flow  relative  to  the  wake --the  wake  points,  in  fact,  convect  with  the 
flow  and  so  the  wake  surface  is  always  aligned  with  the  local  flow.  (The 
entrainment  effect  due  to  turbulent  mixing  is  neglected  for  the  moment.) 

Thus,  Eq.  (5.8)  becomes 

I  I  M  n  .  7  [i]  dS  -  2irMp  +  {|fdS  +  ||M;jii*7[i]dW-0  (5.12) 

S-P  s  w 


In  the  general  case  of  analyzing  the  flow  about  a  given  configuration 
the  doublet  distribution  on  the  surface  is  unknown  while  the  source  distri¬ 
bution  is  determined  directly  by  the  external  Neumann  boundary  condition 
specifying  the  resultant  normal  velocity  at  the  boundary.  The  flow  velocity 
relative  to  the  surface  is, 

V  -  V  -  (5.13) 

where  v  is  the  perturbation  velocity  (Eq.  (5.7)),  and 

-  Vb  >  Q  .  R  -  V.  (5.1A) 

is  the  surface  velocity  relative  to  the  undisturbed  fluid.  V  is  now  mea¬ 
sured  in  an  inertial  frame  (GFF)  which  may  have  an  onset  flow,  v  .  O  is  the 
velocity  of  rotation  of  the  body,  and  R  the  position  of  a  surface  point 
relative  to  the  rotational  axis.  The  normal  component  of  V  is,  from  Eq. 
(5.13), 

V»n-v*n-V  ‘n-V^  (5.15) 

s  N 
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Vjj  is  the  resultant  normal  velocity  at  the  surface.  This  is  usually  zero 
(solid  boundary),  but  it  can  have  a  number  of  nonzero  parts,  e.g.. 


\-vnorm 


(5.16) 


where  VNORM  is  the  user  specified  inflow/outflow  representing  an  engine  in¬ 
let/exhaust  modeling,  and  is  the  boundary  layer  displacement  effect  us¬ 
ing  the  transpiration  technique, 

V-fs<V*) 

where  V  is  the  local  speed  at  the  edge  of  the  boundary  layer  and  S*  is  the 
displacement  thickness.  The  derivative  is  taken  in  the  direction  of  che  lo¬ 
cal  external  flow. 

Using  Eqs.  (5.7),  (5.10)  and  (5.15),  the  source  term  is, 

a-(n*V-+Q*R-vn-n*V  +  VNORM  +  )/4it  (5.18) 

D  A  Dl^ 

When  the  onset  flow  is  nonuniform  in  the  ground- fixed  frame,  the  n  •  term 
for  a  particular  panel  is  time  dependent  as  that  panel  moves  through  space. 
The  local  instantaneous  value  is  picked  up  by  spatial  interpolation  in  the 
defined  onset  flow  according  to  the  location  of  the  panel  center. 

The  wake  development  is  followed  over  a  number  of  time  steps.  At  each 
step  a  new  set  of  wake  points  is  propagated  from  the  wake  shedding  lines, 
each  point  taking  with  it  the  local  doublet  strength  (see  later  under  Kutta 
condition).  Thus,  at  each  time  step  the  wake  doublet  strength  in  Eq.  (5.12) 
is  known.  The  only  unknown  is  the  surface  doublet  distribution  which  is  now 
obtained  as  the  solution  to  Eq.  (5.12). 


Surface  Pressure 


With  the  solution  known,  the  surface  velocities  and  pressures  can  be 
evaluated.  The  tangential  component  of  perturbation  velocity  is  obtained 
from  the  surface  gradient  of  the  potential.  The  normal  component  comes  from 
the  source  term.  The  pressure  coefficient  is  evaluated  using  the  Bernoulli 
equation  for  a  moving  frame. 


C 


P 


V  2 
s 


V*  +  2 

or 


(5.19) 


This  may  be  displayed  in  its  separate  component  parts: 


Dynamic  Term 
Impulsive  Term 
Hydrostatic  Term 


CPD  -  V  2 
s 


V2 


CPI  -  2 


M 


dr 
Fr2 


CPH 
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Here,  the  velocities  are  normalized  by  a  reference  speed,  VREF.  Time  is 
normalized  by  REFL/VREF,  REFL  being  a  reference  length,  e.g.,  propeller  ra¬ 
dius.  f  -  t*VREF/REFL.  The  velocity  potential,  is  normalized  by  VREF 
and  REFL.  The  geometric  height  to  the  free  surface,  z,  is  normalized  by 

REFL,  and  Fr  -  Froude  Number  -  VREF//~gL,  where  L-2*REFL. 


Force  and  Moment 


The  forces  and  moments  are  obtained  by  Integrating  the  pressure  over 
the  surface. 

The  force  coefficient  is, 

‘=F  -  -  H  ”  sSIf  ‘5.20) 

s 


SREF  is  the  reference  area. 

The  moment  coefficient  is. 


Cp  n  A  R  dS 
SREF  L 


(5.21) 


where  L  is  a  reference  length- -usually  the  mean  chord,  CBAR,  for  pitching 
moment  and  SSPAN  (semispan)  for  yawing  and  rolling  moments. 


Kutta  Condition 


If  we  apply  the  condition  of  equal  upper  and  lower  pressure  (using  Eq. 
(5.19))  at  a  thin  trailing  edge,  we  obtain  the  unsteady  Kutta  condition. 


3r 


'W 


+  V, 


3m, 

M  is 


w 


-  0 


(5.22) 


IS 


•u  the  mean  convection  speed  and  the  s  direction  is  the  local  mean  flow 
direction.  My  Is  the  Jump  in  doublet  strength  across  the  trailing  edge,  Eq. 
(5.11),  i.e..  My  Is  the  newly  emerging  wake  strength.  Equation  (5.22)  es¬ 
sentially  states  that  the  rate  of  change  of  circulation  at  the  trailing  edge 
must  match  the  transport  of  circulation  into  the  wake. 
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Wake  Movement 


When  a  solution  has  been  obtained,  velocities  are  computed  at  all  ex¬ 
isting  wake  points  using  a  summation  of  all  singularity  contributions  in  the 
model.  All  wake  points  are  then  propagated  along  the  local  velocity  vector 
for  a  small  time  step.  Simultaneously,  a  new  set  of  wake  panels  is  created 
along  the  separation  lines.  Figure  5.2.  The  current  trailing-edge  doublet 
values  (i.e.,  the  doublet  jump  across  the  wake  at  the  trailing  edge)  is 
transferred  to  the  newly  created  wake  points.  The  doublet  strength  on  each 
wake  point  remains  constant  for  the  remainder  of  the  calculation.  The  wake 
vorticity  effectively  varies  in  time  and  space  according  to  the  local 
stretching  or  contraction  of  the  wake  sheet  as  the  wake  points  convect  at 
the  local  velocities. 

A  new  set  of  doublet  influence  coefficients  is  calculated  for  the  new 
wake  configuration  at  each  time  step. 


64 


USAERO  Configuration  Breakdovm. 


6.0  DISCUSSION  OF  RESULTS 


Five  test  cases  were  chosen  to  validate  different  aspects  of  the  cal¬ 
culation.  These  are  listed  below. 

(1)  Propeller  in  uniform  flow 

(2)  Propeller  in  simple  nonuniform  flow 

(3)  Propeller  in  fully  nonuniform  flow--MPR0P  modeling 

(4)  Unsteady  propeller  flow--USAERO  modeling 

(5)  Simple  propeller-hull  interaction 


6 . 1  Propeller  in  Uniform  Flow 

MPROP  as  modified  and  extended  for  propeller-hull  interaction  was  vali¬ 
dated  against  a  set  of  multi-bladed  propulsors.  Accurate  prediction  of  sur¬ 
face  pressures  is  a  crucial  test  of  the  overall  modeling  capability.  Agree¬ 
ment  with  experiment  for  a  range  of  propeller  models  and  advance  ratios  is 
excellent . 

Figure  6.1  shows  the  paneling  for  a  six-bladed  marine  propeller.  Limi¬ 
ted  blade  surface  pressure  data  was  available  for  this  case  (40) ,  and  the 
comparison  shown  in  Figures  6.2(a)  and  (b)  for  a  typical  radial  station  for 
both  an  on-  and  off-design  condition  shows  good  agreement. 

A  second  propeller  has  been  analyzed  using  MPROP.  This  is  a  flve- 
bladed  propeller  for  which  experimental  results  are  available  at  several 
off-design  conditions  (Hoshino) .  Some  of  these  calculations  are  enclosed  as 
part  of  a  further  demonstration  of  the  code's  capability.  Figure  6.3  shows 
a  five-bladed  propeller.  Pressure  correlation  with  experiment  is  shown  for 
the  advance  coefficients,  J  -  0.3,  0.5  and  0.7,  at  each  of  two  radii  loca¬ 
tions  in  Figures  6.4  through  6.9.  This  propeller  has  a  finite  thickness  at 
the  trailing  edge,  and  the  code  now  has  the  capability  to  treat  correctly 
the  Kutta  condition  for  such  a  geometry.  The  satisfaction  of  the  Kutta  con¬ 
dition  is  achieved  through  iterative  determination  of  the  sectional  circula¬ 
tion.  The  successful  calculations  for  this  propeller  validate  the  newly 
coded  procedure. 

The  paneling  for  a  four-bladed  propeller  is  shown  in  Figure  6.10.  Ex¬ 
tensive  data  on  surface  pressure  is  available  for  this  propeller  tested  in 
Amsterdam  (41).  Figures  6.11  through  6.18  plot  calculated  and  measured  Cp 
at  several  radial  locations  for  the  advance  coefficients,  J  -  0.4  and  J  - 
0.6,  respectively.  Agreement  between  theory  and  experiment  for  this  propel¬ 
ler  is  particularly  good.  The  improved  correlation  is  due  to  the  Y-wake 
model  developed  for  thick  trailing  edges. 
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Figure  6.1.  Blade  Paneling  and  Wake.  Propeller  Model  No.  0123, 
Ship  Research  Institute,  Tokyo,  Japan  (40). 
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Figure  6.3.  Hoshlno  Propeller  Paneling  for  Analysis. 
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Figure  6. A.  Experimental  Correlation,  5-Bladed  Hoshlno  Propeller;  Advance  Ratio 


UPPER  SURFACE 
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Figure  6.5.  Experimental  Correlation,  5-Bladed  Hoshlno  Propeller;  Advance  Ratio  *  0.7,  R  »  10.05. 
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Figure  6.6.  Hoshlno  Propeller  Experimental  Correlation;  Advance  Coefficient  »  0.5,  R  ■  11.80. 
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xperlmenCal  Correlation;  Advance  Coefficient 


UPPER  SURFACE  A  LOUER  SURFACE  X  EXPERinENT 
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Figure  6.8.  Hoshlno  Propeller  Experimental  Correlation;  Advance  Coefficient  -  0.3  (Off  Design) 
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Figure  6.9  .  Hoshlno  Propeller  Experimental  Correlation;  Advance  Coefficient  -  0.3  (Off  Design) 


Figure  6.10.  van  Gent  4-Bladed  Propeller  Paneled  for  Analysis. 
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Experimental  Correlation  for  van  Gent  Propeller  at  Radial  Section,  R  -  10.19,  Advance 
Coefficient  ■  .4. 
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Figure  6.12.  Experimental  Correlation  for  van  Gent  Propeller  at  Radial  Section,  r  -  12.20, 
Advance  Coefficient  -  .A. 
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Experimental  Correlation  for  van.  Gent  Propeller  at  Radial  Section,  r  ■  14.50, 
Advance  Coefficient  ••  .4. 


UPPER  SURFmCE  a  LOUER  SURFACE  X  EXPERIHENT  (41) 
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Experimental  Correlation  for  Van  Gent  Propeller  at  Radial  Section  r  •  19.40, 
Advance  Coefficient  -  .4. 


UPPER  SURFACE  A  LOUER  SURFACE  X  EXPERiriEhT 


Experimental  Correlation  for  van  Gent  Propeller  at  Radial  Section 
Advance  Coefficient  -  0.6. 
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UPPER  SURr'FiCE  A  LOUER  SURFACE  X  EXPERIflEUT  (41) 
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Figure  6.18.  Experimental  Correlation  for  van  Gent  Propeller  at  Radial  Section,  r  ■  19.40 
Advance  Ratio  ••  0.6. 


6 . 2  Propellers  in  Simple  Nonuniform  Flow 


Propellers  operating  in  the  boundary  layer  of  simple  hull  forms  without 
extensive  separation  regions  can  be  very  simply  modeled  with  MPROP.  The 
nonuniformity  can  be  modeled  as  outlined  in  Section  2.2.  The  vorticity  as¬ 
sociated  with  the  shear  flow  is  neglected  in  this  modeling.  The  change  in 
the  boundary  condition  normal  to  the  propeller  surface  is  accurately  model¬ 
ed.  Experimental  results  for  specific  comparison  are  not  available  for  gen¬ 
eral  documentation.  A  typical  calculation  for  a  propulsor  in  radially  non- 
uniform  flow  is  presented.  Figure  6.19  shows  a  front  view  of  a  five-bladed, 
high-skew  propeller  designed  by  Brockett  (P4383).  The  nonuniformity  is  lin¬ 
ear  with  a  lower  free-stream  velocity  (approximately  12%)  near  the  tip  com¬ 
pared  to  the  axial  value.  The  total  pressure  is  assumed  constant  in  the 
radial  direction  and  the  flow  is  nondimensionalized  with  respect  to  the  axi¬ 
al  velocity.  Figures  6.20(a)  and  (b)  show  sectional  Cp  distributions  at  a 
chosen  location  (r  -  4.4)  for  the  basic  propeller  in  uniform  and  nonuniform 
flow,  respectively.  The  change  in  Cp  distribution  due  to  nonuniformity  is 
indicative  of  the  kind  of  detail  that  is  calculated  by  the  code.  The  origi¬ 
nal  propeller  is  believed  to  be  designed  for  a  marginally  nonunifozm  Inflow. 
Analysis  predicts  a  larger  suction  peak  with  nonuniformity.  Further  appli¬ 
cation  of  this  modeling  will  reveal  the  extent  to  which  the  neglect  of  vor¬ 
ticity  is  valid. 


6.3  Propeller  in  Fully  Nonunlform  Flow 

Experimental  results  are  available  for  two  propulsors  in  circumferen¬ 
tially  varying  nonuniform  inflow.  The  nonuniformity  was  simulated  in  the 
test  (40)  to  study  cavitation  under  actual  operating  conditions. 


6.3.1  Frequency  Band  Analysis 

Both  propellers  have  been  analyzed  with  the  frequency  band  theory  out¬ 
lined  in  Section  5.2.  In  Figure  6.21  the  experimental  wake  fractions  are 
reproduced  from  Ref.  40.  The  nonuniformity  was  represented  by  a  Fourier  co¬ 
sine  series.  The  second  harmonics  or  the  cos  (2wt)  term  is  the  most  domi¬ 
nant.  Figure  6.22  shows  the  comparison  between  calculated  and  measured 
pressure  at  a  radial  section,  r  -  0.77R  for  J  -  .9.  The  experimental  pres¬ 
sure  data  is  available  at  a  gauge  located  at  20%  downstream  of  the  leading 
edge.  The  dotted  line  in  Figure  6.22  shows  a  three -component  solution  for 
the  velocity  potential.  The  calculated  pressure  is  in  good  agreement  with 
experiment.  The  slight  overprediction  of  the  peak  pressure  is  not  fully 
understood.  It  is  believed  that  the  theory  forces  a  continuous  analytical 
function  for  the  velocity  potential  which  introduces  inaccuracies  in  the 
time  derivative  of  the  potential  which  influences  the  pressure  calculated. 
A  five-term  solution  is  also  represented  in  Figure  6.22  by  the  solid  line. 
The  change  in  solution  is  not  significant  because  the  amplitude  of  the  high¬ 
er  harmonics  is  very  small.  Full  details  of  the  experimental  setup  are  not 
available.  There  is  a  change  in  hydrostatic  pressure  across  the  radius  of 
the  propeller  which  has  not  been  accounted  for  in  the  theory. 
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Figure  6.19,  Propeller  with  a  Skew  of  72  Degrees  In  Uniform  Flow. 
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Reference  40 


EXPERIMENT  (40)  -  3  TERMS  1  CALCULATION 

-  5  TERMS  J  MPROP 
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Figure  6.22.  Pressure  Correlation  on  a  Slx-Bladed  Propeller  in  Nonunlfo 
x/c  ■■  0.2,  Suction  Side. 


The  nonuni formlty  associated  with  a  five-bladed  Hoshino  propeller  is 
shown  in  Figure  6.23.  The  description  of  nonuniformity  is  less  precise  than 
the  previous  case.  The  contour  plots  were  ditigized  and  represented  by  an 
8 -term  Fourier  series.  The  wake  fraction  locally  drops  to  0.7  in  a  narrow 
region  around  the  vertical  axis.  This  makes  it  necessary  to  have  many  terms 
before  an  acceptable  solution  can  be  calculated.  Figure  6.24  shows  the  cal¬ 
culated  results  plotted  against  experimental  results  taken  from  Ref.  42. 
Experimental  probes  were  located  at  a  radial  section  0.75  R  outboard.  The 
suction  side  probe  was  located  at  10%  downstream  of  the  leading  edge  while 
the  pressure  side  probe  was  located  40%  downstream  of  the  leading  edge.  The 
calculated  pressures  are  in  good  agreement  for  the  complex  inflow.  The  ba¬ 
sic  assumptions  and  approximations  of  the  theory  may  be  considered  indirect¬ 
ly  validated  by  these  calculations, 

6.3.2  Time -Dependent  Analysis 

The  existing  time -dependent  program,  USAERO,  described  in  Section  5.3, 
was  applied  to  the  SRI  MP0123  propeller  (40) . 

6. 3. 2.1  Uniform  Flow 

First,  as  in  the  case  of  the  frequency  band  analysis,  the  uniform  flow 
case  (40)  was  examined  to  establish  a  basic  level  of  accuracy.  A  relatively 
low  panel  density  representation  of  the  blades  was  used  (Figure  6.25),  i.e,, 
208  panels  on  each  blade  in  a  16  (chordwise)  x  12  (radially)  array  plus  a  2 
X  8  array  on  the  tip.  A  simple  hub  was  also  included  to  improve  the  blade 
root  representation.  The  propeller  was  started  impulsively  from  rest  and 
run  for  three  revolutions  at  10  rps  with  12  steps  in  each  of  the  first  two 
revolutions  and  24  steps  in  the  third.  These  runs  were  made  with  different 
onset  flows  to  give  advance  ratios  of  0.7,  0.9  and  1.1. 

Figure  6.26(a)  shows  a  general  view  of  the  pressure  coefficient  distri¬ 
bution  after  three  revolutions  at  J  -  0.7.  (Note  that  hydrostatic  pressure 
is  included  and  therefore  the  Cp  varies  with  azimuthal  location.)  Stream- 
wise  lines  indicate  the  computed  location  of  the  wake  surfaces.  A  buttline 
cut  through  the  wake  indicates  a  higher  rate  of  convection  in  the  root  re¬ 
gion  compared  with  the  tip  region.  Figure  6.26(b).  The  starting  vortex  from 
the  impulsive  start  has  convected  about  2  1/2  diameters  downstream.  A  sta¬ 
tion  cut  through  the  wake  at  one  diameter  downstream  from  the  propeller. 
Figure  6.26(c),  shows  a  large  middle  region  in  essentially  "solid-body" 
swirl  and  a  thin  outer  region  where  the  tip  vortex  flow  is  becoming  sheared 
due  to  lower  swirl  rate  there.  An  inner  region  involved  with  the  hub  root 
vortex  is  also  evident. 

Figure  6.26(d)  shows  the  time  history  of  x-wise  force  (-CD)  and  normal 
(CMX)  coefficients  from  the  USAERO  calculation.  They  converge  to  steady 
state  values  in  12  steps,  i.e.,  one  revolution.  The  wake  has  convected  by 
less  than  one  diameter  in  that  time.  Reduced  to  the  standard  thrust  (KT) 
and  torque  (KQ)  coefficients,  these  and  similar  values  obtained  at  J  -  .9 
and  1.1  compare  very  favorably  with  the  experimental  measurements,  Figure 
6.26(e),  in  view  of  the  relatively  low  panel  density  used.  Viscous  effects 
were  not  included  in  the  calculation  and  so  the  slightly  higher  thrust  value 
obtained  is  to  be  expected.  However,  the  slightly  higher  torque  coefficient 
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Figure  6.23.  Wake  Fraction  Contours  for  Nonuniform  Inflow,  Reproduced 
from  Reference  kl. 
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omparlson  between  8-Term  Solution  and  Experiment  (Ref.  42). 
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6  25  L/SAKKO  Paneling  for  the 


Figure  6.26.  USAERO  Calculations  for  the  SRI  MP0123  6-Bladed  Propeller  in  Unifor 


MP0123 


Figure  6.26.  Continued. 


obtained  from  the  calculation  is  not  consistent  with  the  expected  correction 
due  to  viscous  effects;  it  is  likely  that  the  calculation  has  underestimated 
the  leading-edge  thrust  because  of  poor  panel  resolution  around  the  leading 
edge.  A  higher  panel  density,  at  least  in  the  chordwise  direction,  would  be 
essential  in  order  to  predict  torque  correctly.  This  would  also  be  a  re¬ 
quirement  when  including  boundary  layer  calculations  in  the  analysis. 

Figures  6.26(f),  (g)  and  (h)  show  a  very  reasonable  correlation  between 
the  calculated  pressure  distribution  and  measured  values  at  .77  radius  for 
advance  rations  of  .7,  .9  and  1.1,  respectively.  The  normalizing  velocity 
for  the  experimental  pressure  coefficient  is  the  local  inflow  velocity  plus 
an  induced  term  computed  by  a  lifting  surface  theory.  The  latter  term  is 
unknown  here  but  shjuld  be  small  at  the  advance  ratios  considered.  For  the 
present  calculation,  therefore,  the  normalizing  velocity  for  pressure 
coefficient  is  just  the  local  inflow  velocity,  i.e.,  the  combination  of  the 
onset  flow  and  the  velocity  due  to  rotation  at  .77  radius. 

The  results  discussed  above  establish  a  basic  level  of  accuracy  for  the 
code.  Although  the  application  of  a  time-stepping  method  to  the  analysis  of 
propellers  in  uniform  flow  would  normally  be  regarded  as  "overkill",  it  is 
interesting  to  nota  that  for  the  advance  ratios  covered,  steady  state  solu¬ 
tions  were  reached  in  less  than  15  minutes  on  an  SGI  IRIS  Power  workstation. 
Moreover,  the  calculations  include  the  natural  convection  of  the  wakes. 

6. 3.2.2  Nonuniform  Flow  Analysis 

In  the  nonuniform  onset  flow  case  the  experimental  (40)  wake  fractions 
(Figure  6.21)  were  digitized  and  placed  directly  into  the  USAERO  flow  cor¬ 
responding  to  J  -  .9.  The  resulting  axial  onset  flow  is  displayed  in  a 
spectrum  plot  in  Figure  6.27.  Red  indicates  the  low  flow  region  while  blue 
represents  the  undisturbed  flow.  The  velocity  is  normalized  by  the  speed  of 
rotation  at  .77  of  the  propeller  radius.  The  circular  survey  plane  has  a 
diameter  1.2  times  that  of  the  propeller  and  is  placed  just  upstream  of  the 
propeller . 

Figure  6.28  shows  a  general  view  of  the  propeller  with  a  moderate  panel 
density;  480  panels  were  used  on  each  blade  in  a  30  x  15  array  plus  2  x  15 
on  the  tip.  This  is  approximately  double  the  panel  density  used  for  the 
uniform  flow  case  but  is  still  about  half  the  density  used  earlier  in  the 
MPROP  frequency  band  analysis.  The  propeller  was  started  impulsively  from 
rest  to  10.8  rps  at  J  -  .9  and  was  run  for  three  revolutions  with  24  steps 
per  revolution.  Only  results  from  the  last  complete  revolution  were  passed 
to  the  plot  file.  Figure  6.28  shows  the  computed  surface  pressure  distribu¬ 
tion  at  the  end  of  the  third  revolution.  The  pressure  values  include  the 
dynamic,  impulsive  and  hydrostatic  terms.  Note  the  relatively  high  leading- 
edge  suction  levels  (red)  on  the  blades  near  the  top  and  bottom  azimuthal 
positions  compared  with  lower  suction  levels  (green)  on  the  two  blades  in 
the  high  flow  region:  the  rotational  velocity  combined  with  the  full  onset 
flow  produces  a  lower  local  angle  of  attack  than  when  combined  with  the 
smaller  onset  flow  in  the  "wake".  The  propeller  wake  itself  is  displayed  by 
streamwise  lines.  These  indicate  the  computed  locations  of  the  wake  doublet 
sheets  which  are  propagated  from  each  blade  trailing  edge  and  convected 
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Figure  6.27.  Resulting  Onset  Flow  Computed  on  a  1.2  D  Disc  Just  Upstream  of  Propeller 
(Vx  Normalized  by  Speed  of  Rotation  at  .77  R) . 


General  View  of  the  MP0123  Propeller  in  Nonuniform  Flow  Showing  the 
USAERO  Computed  Surface  Pressure  Distribution. 


downstream  step-by-step  at  the  local  velocities.  In  this  case  the  convec¬ 
tion  rate  is  slower  at  the  top  and  bottom  azimuthal  stations  compared  with 
those  at  the  sides  and  so  there  is  a  tendency  to  bunch  up  near  the  vertical 
plane.  Figures  6.29(a)  and  (b)  show  two  waterline  cuts  at  Z  -  3.0  and  7.0 
(.24  and  .56  of  blade  radius,  respectively)  in  the  low  flow  region  and  indi¬ 
cate  a  smooth  wake  propagation  from  the  blade  trailing  edges  at  the  various 
conditions  captured  by  the  cuts  on  different  blades. 

Because  of  earlier  problems  encountered  in  MPROP  in  regard  to  the  Kutta 
condition  in  the  blade  root  region,  the  computed  pressure  distribution  was 
also  inspected  at  the  Z  -  3.0  waterline  cut  (Z/R  -  .24;  the  hub  is  at  .18). 
Figure  6.30  indicates  a  very  good  behavior  of  the  unsteady  Kutta  condition 
employed  in  USAERO.  The  impulsive  term  (included  in  the  total)  is  also  dis¬ 
played  (at  the  same  scale)  to  indicate  the  magnitude  of  the  unsteady  effect 
at  this  position. 

Typical  chordwise  pressure  distributions  at  approximately  .77  radius 
are  given  in  Figure  6.31(a),  (b) ,  (c)  and  (d)  for  a  blade  at  azimuthal  sta¬ 
tions  0,  90,  180  and  270  degrees,  respectively.  They  are  all  plotted  to  the 
same  Cp  scale  to  Indicate  visually  the  range  of  pressure  variation  in  one 
revolution.  They  all  indicate  a  good  behavior  of  the  unsteady  Kutta  condi¬ 
tion  over  a  wide  range  of  conditions.  The  impulsive  term  is  included  as 
well  as  the  complete  pressure  coefficient.  A  load  reversal  is  evident  at 
the  90*  and  270*  positions,  i.e.,  in  the  full  onset  flow.  This  is  consis¬ 
tent  with  the  general  view  in  Figure  6.28*. 

Experimental  pressure  measurements  (40)  were  taken  at  five  chordwise 
positions  at  .77  radius  on  the  blade.  Positions  A,  B,  C  and  D  are  on  the 
"upper"  or  "suction"  side  at  .2,  .4.  .6  and  .8  of  the  chord,  respectively, 
from  the  local  leading  edge  and  position  E  is  at  .4  of  the  chord  on  the 
"pressure"  side.  Comparisons  between  measured  and  calculated  pressures  for 
one  complete  revolution  are  shown  in  Figure  6.32(a),  (b),  (c) ,  (d)  and  (e) 
for  positions  A,  B,  C,  D  and  E,  respectively.  The  calculations  include 
curves  for  a  lower  panel  density  case  which  used  240  panels  in  a  16  x  14  ar¬ 
ray  plus  2  x  8  on  the  tip.  At  position  A  (-2c)  there  is  very  little  effect 
of  panel  density  but  the  calculated  pressures  are  generally  more  negative 
than  the  experimental  values.  (The  frequency  band  analysis  tended  to  be 
more  positive,  Figure  6.22.)  Although  the  main  features  of  the  pressure 
variation  are  captured  by  the  calculation  there  is  a  general  appearance  of  a 
phase  lead.  Basically,  the  calculation  follows  the  experiment  fairly  close¬ 
ly  in  the  region  of  decreasing  peak  suction  (i.e.,  when  entering  the  high 
onset  flow  region),  but  it  leads  the  experiment  on  the  increasing  peak  suc¬ 
tion  region  (entering  the  low  onset  flow  region).  The  effect  of  viscosity 
could  be  a  factor  here;  boundary  layer  effects  were  not  included  in  the  pre¬ 
sent  calculations.  At  position  B  (-4c),  the  calculation  is  again  more  nega¬ 
tive  chan  the  measured  values  but  in  this  case  the  discrepancy  is  ar.  almost 
constant  10  mm.  of  Hg.  In  this  case  the  main  features  of  the  pressure  vari¬ 
ation-  -which  is  less  extreme  than  at  postlon  A- -are  captured  by  the  calcula¬ 
tion  and  are  more  or  less  in  phase.  The  lower  panel  density  calculation  at 
this  position  is  again  fairly  close  to  the  high  density  case;  the  main  devi¬ 
ations  occur  near  the  bottom  azimuthal  position  where  the  blade  is  leaving 
the  low  onset  flow  region.  This  deviation  of  the  low  panel  density  case  is 
somewhat  more  obvious  for  position  C  (-6c)  at  the  same  azimuthal  station. 
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WATERLINE  CLIT  SRI  MP0123  6-BLADED  PROPELLER  in  non-uniform  flow  J=0.9 
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PRESSURE  HISTORY  COMPARISON  (position  A) 
exp.  versus  cole,  (high  ond  low  ponel  density) 


PRESSURE  HISTORY  COMPARISON  (position  B) 
exp.  versus  cole,  (high  ond  low  ponel  density) 
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PRESSURE  HISTORY  COMPARISON  (position  C) 
exp.  versus  cole,  (high  ond  low  ponel  density) 
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PRESSURE  HISTORY  COMPARISON  (position  E) 
exp.  versus  cole,  (high  ond  low  ponel  density) 
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azimuth  Idegrees) 

SRI  MP0123  PROPELLER,  USAERO  CALCULATION 


Again,  the  high  density  case  follows  the  experimental  pressure  variation  at 
position  C  very  closely  but  with  a  general  negative  offset  of  about  10  mm. 
Hg. 


At  position  D,  Figure  6.32(d),  the  high  density  case  again  follows  the 
experimental  pressure  variation  very  closely  but  with  a  small  negative  off¬ 
set.  The  variation  of  pressure  here  is  relatively  benign  compared  with  that 
at  position  A.  The  shift  of  the  low  panel  density  calculation  when  leaving 
the  low-onset-flow  region  is  even  more  pronounced  than  at  positions  B  and  C. 
This  is  caused  by  a  small  shift  in  the  phase  of  the  impulsive  term  in  the 
lower  density  case  compared  with  that  of  the  higher  density  paneling. 

Finally,  at  position  E  (-4c  on  the  "pressure"  side)  the  calculation 
again  captures  the  main  features  of  the  measured  pressure  variation.  This 
time  there  is  a  tendency  for  a  positive  shift  in  pressure.  There  is  little 
effect  of  panel  density  evident  at  this  position.  Overall,  the  tendency  for 
a  negative  pressure  shift  on  the  upper  surface  and  positive  shift  on  the 
lower  surface  are  consistent  with  viscous  effects  and  so  it  would  be  expect¬ 
ed  that  the  inclusion  of  boundary  layer  calculations  in  the  unsteady  calcu¬ 
lation  would  improve  the  correlation  with  experiment. 

It  is  instructive  to  examine  the  detailed  breakdown  of  the  pressure 
variation  into  its  component  parts,  i.e.,  dynamic,  impulsive  and  hydrostat¬ 
ic,  as  outlined  in  Eq.  (5.19).  Figure  6.33(a)  shows  the  details  for  posi¬ 
tion  A.  The  dynamic  term  has  a  peak- to- trough  variation  of  about  70  mm.  Hg 
and  lags  behind  the  impulsive  term  by  about  75*  of  azimuth.  The  impulsive 
term  has  a  peak- to -trough  variation  of  about  40  mm.  Hg.  Both  the  dynamic 
and  impulsive  terms  have  two  cycles  per  revolution  while  the  hydrostatic 
term,  which  has  a  peak- to -trough  variation  of  about  16  mm.  Hg,  has  one 
cycle.  With  such  large  variations  in  the  major  terms  and  in  view  of  the 
phase  difference  between  them,  it  is  remarkable  that  the  combined  total  is 
as  close  to  the  experiment  as  it  is  (Figure  6.32(a)).  It  is  even  more  re¬ 
markable  at  position  D  (Figure  6.33(b))  where  the  resultant  pressure  varia¬ 
tion  is  relatively  benign.  Here  the  peak- to -trough  variations  in  the  dynam¬ 
ic  and  impulsive  terms  are  50  and  70  mm.  Hg,  respectively,  and  they  are  es¬ 
sentially  90*  out  of  phase. 

Figure  6.34  displays  the  dynamic  and  impulsive  pressure  variations  in 
relation  to  the  onset  flow  sensed  at  position  A.  If  we  think  in  terms  of 
the  angle-of-attack  variation  due  to  the  nonuniform  onset  flow,  i.e.,  a 
small  in  the  high  flow  region  and  a  large  in  the  low  flow  region,  then  the 
dynamic  pressure  is  lagging  a  by  about  25*  while  the  impulsive  term  is  lead¬ 
ing  by  about  50* . 

While  it  is  true  that  the  phasing  between  the  dynamic  and  impulsive 
pressure  terms  is  strongly  connected  to  the  reduced  frequency  of  the  effec¬ 
tive  a  variation  at  each  station,  it  is  not  clear  to  what  extent  the  fea¬ 
tures  of  blade  design,  i.e.,  rake,  skew,  twist,  camber  and  thickness  distri¬ 
bution,  affect  the  relationship.  The  time-stepping  approach  provides  a 
means  of  examining  these  factors  with  a  view  towards  achieving  a  blade  de¬ 
sign  which  minimizes  the  load  variations  for  a  given  performance  on  a  given 
hull  (i.e.,  in  the  presence  of  a  given  nonuniform  onset  flow). 
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PRESSURE  HISTORY  DETAILS  (position  A) 
showing  dynomic,impulsive  ond  hydrostatic  terms 


PRESSURE  HISTORY  DETAILS  (position  0) 
showing  dynomic, impulsive  ond  hydrostotic  terms 


PRESSURE  PHASING  WITH  ONSET  FLOW  (ot  station  A) 

showing  dynomic  ond  impulsive  terms 
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6.4  Simple  Propeller-Hull  Interaction 

Figure  6.35  shows  a  four-bladed  propulsor  operating  in  the  wake  of  a 
long  hull.  The  hull  geometry  is  similar  to  that  described  in  Ref.  33  as 
Afterbody  1.  The  propeller  geometry  is  derived  from  Ref.  41  with  appropri¬ 
ate  scaling  to  match  the  maximum  hull  diameter.  The  calculation  of  viscous 
flow  on  the  hull  for  a  Reynolds  number  of  4.56  million  was  accomplished  with 
RANSTERN  as  described  in  section  4.2.  The  calculation  of  propeller  flow 
with  MPROP  in  the  absence  of  the  hull  has  been  extensively  correlated  with 
experiment  as  discussed  in  Section  6.1.  In  this  section  the  effect  of  hull 
wake  flow  on  the  propeller  is  considered. 

Figure  6.36  shows  the  axial  velocity  distribution  as  a  function  of 
radius.  "Full"  velocity  profile  refers  to  the  case  where  the  propeller  di¬ 
ameter  is  equal  to  the  maximum  hull  diameter.  "Partial"  velocity  profile 
refers  to  the  case  where  the  hull  diameter  is  a  fifth  of  the  propeller  di¬ 
ameter.  The  nonuniformity  of  the  flow  is  limited  to  the  hub  region  in  the 
latter  case  with  uniform  flow  over  most  of  the  propeller.  The  hull  stern 
velocity  profiles  were  calculated  using  RANSTERN  as  previously  described. 
Pressure  contours  for  two  advanced  coefficients  (J  -  0.4,  and  0.5)  were 
calculated  with  MPROP  for  prescribed  nonuniformity.  Figures  6.37(a)  through 
(J)  present  spectrum  plots  of  calculated  pressures  for  a  variety  of  cases 
ranging  from  uniform  flow  (no  hull  interaction)  to  full  presence  of  hull 
wake.  For  uniformity  the  scale  chosen  for  all  pressure  contours  is  -4.0  to 
4.0.  For  this  choice  of  pressure  range,  the  contours  on  the  "pressure"  side 
of  the  propeller  (Figures  6.37(b),  (d),  (f),  (h)  and  (j))  do  not  show  much 
variation.  A  comparison  of  Figure  6.37(a)  with  Figure  6.37(c)  indicates  the 
effect  of  hull  wake  on  the  calculated  pressure  distribution  for  the  J  -  0.4 
case.  A  similar  comparison  of  Figure  6.37(g)  with  6.37(i)  indicates  the  in¬ 
creased  effect  of  hull  wake  on  the  propeller  for  higher  advance  coeffici¬ 
ents.  Figures  6.37(e)  and  (f)  show  a  case  with  very  slender  hull,  i.e., 
partial  nonuniformity.  The  changes  to  the  pressure  levels  as  compared  with 
the  uniform  flow  case  would  be  expected  to  be  small  for  this  case.  This  ex¬ 
pectation  is  borne  out  in  general.  However,  a  particular  feature  of  the  ef¬ 
fect  of  the  reduced  hull  flow  on  the  propeller  in  the  region  near  the  hub 
actually  shows  a  higher  circulation  and  greater  suction  pressure.  This  is 
due  to  the  increased  angle  of  attack  of  the  flow  relative  to  the  local  suc¬ 
tion. 
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Figure  6.37.  Propeller  Pressure  Contours  Calculated  by  MPROP  to  Highlight  Viscous  Hull  Interactions. 


Figure  6.37.  Continued. 
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Figure  6.37.  Continued. 
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Figure  6.37.  Continued. 
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(e)  Slender  Hull;  J  =  0.4,  Suction  Side 


Figure  6.37.  Continued. 
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Figure  6.37.  Ccntinued. 


Figure  6.37.  Continued. 


Figure  6.37.  Continued. 


135 


Figure  6.37.  Concluded. 


7.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  complex  problem  of  propeller/hull  interaction  has  been  examined  in 
separate  parts,  i.e.,  (i)  a  propeller  operating  in  a  flow  generated  by  a 
hull,  and  f^ii)  a  hull  operating  in  the  effect  of  a  propeller. 

Details  of  velocity  profiles  in  the  hull-generated  viscous  flow  ap¬ 
proaching  the  propeller  location  were  computed  using  the  Reynolds -Averaged, 
Navier-Stokes  program,  RANSTERN. 

The  basic  flow  characteristics  of  a  hull  are  treated  adequately  using 
VSAERO,  a  panel  method  with  coupled  boundary  layer  calculations.  Example 
calculations  of  a  propller  operating  in  such  a  nonuniform  flow  were  demon¬ 
strated  for  the  axisymmetric  case  using  the  marine  propeller  version  of 
VSAERO,  MPROP,  with  plausible  results.  The  code  was  also  used  to  demon¬ 
strate  the  simple  (axisymmetric)  case  of  a  hull  operating  in  the  presence  of 
a  propeller. 

The  MPROP  program  was  also  used  to  explore  calculations  of  radially 
varying  onset  flows  for  a  propeller/hub/shroud/stator  configuration.  Again, 
plausible  results  were  obtained  but  comparisons  with  measurements  are  needed 
to  complete  the  study. 

The  more  general  condition  of  a  propeller  operating  in  a  nonaxisymmet- 
ric  hull  wake  is  more  challenging  because  the  propeller  sees  time -dependent 
conditions  as  it  rotates  in  the  arbitrary  inflow  distribution.  Two  ap¬ 
proaches  were  applied  to  this  problem,  both  with  very  encouraging  results. 
The  first  approach  used  a  frequency-band  analysis  of  the  azimuthal  variation 
of  the  inflow  and  solved  for  doublet  Fourier  coefficients  in  the  MPROP  pro¬ 
gram.  This  very  cost  effective  analysis  gave  results  which  compare  very 
well  with  experimental  measurements  even  though  certain  wake  assumptions  are 
needed  in  the  model.  A  second,  more  time  consuming,  but  still  practical 
calculation  was  conducted  with  USAERO,  a  time- stepping  panel  method,  which 
rotates  a  fully  paneled  propeller  in  finite  steps  through  the  nonuniform 
flow.  The  output  from  this  is  a  detailed  time  history  of  velocities,  pres¬ 
sures  and  integrated  forces  and  moments.  This  approach  also  gave  good 
agreement  with  experimental  measurements;  moreover,  it  provided  a  detailed 
breakdown  of  the  pressure  components,  i.e.,  the  dynamic,  impulsive,  and 
hydrostatic  terms,  plus  their  varying  phase  relationships  for  different 
positions  across  a  blade  chord.  Unlike  the  frequency-band  analysis,  the 
time -stepping  program  requires  no  assumptions  as  to  wake  position  or  peri¬ 
odic  strength  variation  since  the  wake  is  generated  step  by  step  according 
to  local  conditions. 


136 


7 . 1  Recominendations 


(1)  In  view  of  its  capability  and  cost  effectiveness,  the  MPROP  version  of 
VSAERO  should  be  developed  further  to  improve  its  usability  and  applic* 
ability.  Certain  refinements  in  the  wake  model  would  improve  its 
treatment  of  propellers  at  low  advance  coefficients.  Also,  the  model¬ 
ing  of  the  periodic  wake  strength  in  the  frequency  band  analysis  of 
nonuniform  flows  could  be  improved.  The  basic  program  should  also  be 
updated  to  take  advantage  of  more  recent  developments  in  the  VSAERO 
program,  especially  in  regard  to  the  coupled  boundary  layer  analysis. 
Overall,  more  comparisons  with  measured  data  are  needed  in  order  to  es¬ 
tablish  its  range  of  validity. 

(2)  A  hybrid  coupling  between  RANSTERN  and  VSAERO  should  be  pursued  with  a 
view  to  providing  a  cost  effective  zonal  approach  to  computing  viscous 
velocity  profiles  in  the  propeller  region. 

(3)  The  USAERO  program  should  be  applied  to  more  cases  for  comparison  with 
measured  data.  In  particular,  boundary  layer  calculations  should  be 
included.  The  boundary  layer  routines,  although  coupled  in  an  earlier 
version  of  USAERO,  have  been  temporarily  out  of  action  during  an  exten¬ 
sive  rewrite  and  restructuring  of  the  basic  code.  These  routines  are 
about  to  be  reactiveated. 

The  very  general  modeling  capability  of  the  program  offers  a  potential¬ 
ly  broad  range  of  propeller/hull  interaction  problems  that  could  be 
employed:  for  example, 

(a)  a  propeller  in  general  orientation/motion,  e.g.,  in  yaw,  in  surge, 
heave  and  sway,  etc.,  simulating  the  effect  of  hull  motions  on  the 
rotation  propeller; 

(b)  a  propeller  with  hub,  shroud,  and  stators; 

(c)  complete  propeller/hull  configuration  for  simultaneous  solution  of 
the  mutual  interaction  effects  including  periodic  pressures  on  the 
hull  surface;  and 

(d)  a  propeller  operating  near  a  free  surface. 


USAERO  should  also  be  used  to  study  in  more  detail  the  complex  phasing 
of  the  dynamic  and  impulsive  pressure  terms  on  a  propeller  operating  in  a 
nonunifom  flow.  The  objective  would  be  to  Improve  the  understanding  of  the 
influence  of  various  blade  parameters  (skew,  rake,  twist,  camber,  etc.)  on 
this  phasing  in  order  to  reduce  the  resultant  load  oscillations  for  a  given 
performance  in  blade  design. 
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